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FOREWORD 

This document is Volume I of the final report for Contract NAS8-31310. 
Volume II is printed under separate cover and is entitled "Solid Rocket Booster 
Thermal Radiation Model — Volume II — User’s Manual," LMSC-HREC TR 
D496763-II. This report was prepared by personnel of the Thermal &; Fluid 
Physics Group, Engineering Sciences Section, Lockheed-Huntsville Research 
& Engineering Center. The contract period of performance was from 20 
January 1975 through 20 March 1976. The work was administered under the 
technical direction of Mr. William C. Claunch of the Structures and Propulsion 
Laboratory, NASA -Marshall Space Flight Center. 
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Symbol 

English 

A 

A/E 



e 


F 

F(x) 

H. 

1 

h 

I 

J 


K 


K 


a 


N 

N 

e 



n 


1,2,3 


NOMENCLATURE 

Definition 


area 

absorption to extinction ratio 
quadric coefficients 
param eter 

black body emissive power 
energy per bundle 
view factor 

probability density function 
elevation of axial region boundaries 
elevation along plume axis 

radiation intensity (also axial region boundary index) 

volumetric emission of radiant energy (also radial 
region boundary index) 

particle -wavelength size parameter 

linear absorption coefficient 

length 

particle number density 
sample size (number of emissions) 
number of hits on a surface 
direction cosines of line of sight 
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P 

P 


1. 2. 3 


P 

c 

P. 

1 

P 

00 

P(x) 

^a, e, scat 

q 

R 


R 


ex 


r 


P 

S 

S 

P 

S(t},6) 

T 

U 


'^1, 2, 3 
V 

^1, 2. 3 
Z 


Greek 

p 

y 


points on a surface 

pressure or property 

chamber pressure 

intersection point 

ambient pressure 

probability distribution function 

absorption, extinction, scattering efficiencies 

heating rate 

distance between plume axis and line of sight in plane 
at right angles to plume*s axis passing through cone 
vertex 

nozzle exit radius 

radial coordinate (cylindrical coordinate) 

particle radius 

length along line of sight 

probable path length 

scatter function 

temperature 

random number 

direction cosines of surface unit normal 
volume 

principal coordinates (Cartesian coordinates) 

distance along plume axis from base plane (cylindrical 
coordinate ) 

angle between r^ and r^ 

cone half angle for radial region boundaries 
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n 

9 

4 > 

4 ^ 

a 

a, e, scat 

T 

0) 

X 

Subscripts 

a 

b 

g 

P 

scat 

X 


polar angle of line of sight (spherical coordinate) 
azimuth angle of line of sight (spherical coordinate) 
azimuth angle of event site (cylindrical coordinate) 
gimbal angle in yaw plane 

gimbal angle in pitch plane (also Stefan -Boltzmann 
con stant) 

absorption, extinction, scattering cross section for 
AjC^O^ particles 

optical depth 

solid angle 

wavelength 

absorbtion 

blackbody 

gas 

particle 

scattering 

wavelength 
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1. INTRODUCTION 

Thermal radiation from exhaust plumes of solid rocket boosters (SRB) 
on the Space Shuttle is a major contributor to the heat load incident on the 
base region of the Space Shuttle vehicle. The magnitude of this heating rate 
is expected to be comparable to that caused by the Space Shuttle main engines 
(SSME). Accurate estimates of SRB thermal radiation is necessary to opti- 
mally design the thermal protection system (TPS) for the ascent phase of the 
Space Shuttle trajectory. These estimates often rely on analytical models of 
the thermal radiation emitted by plumes. 

SRB plumes are composed of a dispersion of particles distributed 

throughout the gaseous products of combustion of the solid propellant. 
particles are present in SRB plumes because alumina is added to the solid 
propellant to increase the specific impulse of the booster. The particles may 
exist in either the solid or liquid phase. In the solid phase, P^^^ticles 

are almost pure scatterers of electromagnetic radiation. However, in the 
liquid phase Alt^O^ has a relatively large absorbtion cross section and con- 
tributes significantly to the thermal radiation of the SRB plume. Certain 
gaseous constituents present in SRB plumes are also strong emitters of 
radiant energy. So, a realistic model for thermal radiation from SRB plumes 
must consider coupled radiation from particles as well as emission 

from the principal emitting gases, namely CO, CO^, and HCf. 

At the present time, much of the thermal environment calculations in 
the base region of the Space Shuttle which account for radiant heating due to 
SRB plumes are based on a view factor calculation to surfaces in the base 
region from a diffuse conical surface surrounding the plume. These view 
factor calculations are used in conjunction with an assumed plume emissive 
power in order to calculate radiant heating rates. Assumptions are also 
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necessary regarding an appropriate value for the half angle to be used for 
the diffuse conical surface representing the SRB plume. Unfortunately, view 
factor methods for calculating thermal environments due to SRB plumes 
cannot account for influences on plume radiant heating of such parameters 
as: 

• Chamber pressure 

• Propellant composition 

• Altitude 

• Approach Mach number 

• Afterburning, and 

• Spatial variations in plume properties, including: 

(1) gas and particle temperature; (2) gas 

pressure; (3) Ai 203 particle number density; (4) 
gas mole fractions; and (5) scatter and absorption 
cross sections. 

This report documents the development and description of a SRB plume 
thermal radiation model which accounts for the influence of the above param- 
eters on the magnitude and distribution of thermal radiation emanating from 
SRB plumes. The Monte Carlo method is used to solve the radiation transport 
problem to provide more flexibility in modeling plume structure. Property 
definition in the plume may be provided by flowfield calculations or, if avail- 
able, from experimental data. Spatial distributions of plume properties are 
approximated by arranging homogeneous regions throughout the plume (in 
which plume properties are defined) in a manner consistent with plume prop- 
erty variations. 

The Space Shuttle will have two SRBs astride the vehicle that will con- 
tribute to the base heating. These plumes may intersect close to the base 
region at high altitudes and may also gimbal one independent of the other. 
Plume intersection and plume gimbal angle can have a considerable effect 
on local heating rates. In view of this expected character of SRB thermal 
radiation, the SRB radiation model described herein considers dual SRB 
plumes. In this manner, the attenuation of radiant energy from one plume 


1-2 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D496763-I 

on its neighboring plume (or shadowing effect) is accounted for in the model's 
calculations of the distribution of thermal energy leaving the plume. 

Space Shuttle surface geometries are modeled as combinations of 
quadric surfaces. Specific surface geometries that can be considered are; 

• Planar (parallelegram or annular) 

• Cones and frustums 

• Cylinders 

• Ellipsoids and spheres, and 

• Paraboloids. 

These quadric surfaces can be defined in a Cartesian coordinate system by 
means of a single polynomial expression. Models of specific Space Shuttle 
surface geometries can be easily constructed with combinations of these 
surfaces. The individual surfaces can be subdivided to provide more dis- 
crete definition of local geometries. Each subdivision of every surface has 
a unique mathematical definition. This surface characteristic is used in 
conjunction with the Monte Carlo solution to the radiation transport problem 
so that radiation heating rates can be calculated to each discrete subsurface. 
The Monte Carlo method also enables easy and accurate predictions of the 
effects of shading of one surface by another on the magnitude of radiant heat 
loads. 

The SRB thermal radiation computer code also provides the capability 
to calculate view factors between one selected surface on the Space Shuttle 
and the SRB plumes as well as all other surfaces of interest on the Space 
Shuttle. This capability will allow detailed analysis of SRB thermal radia- 
tion at extremely localized areas of the Space Shuttle where the effects of 
surface shading are considerable. 

The SRB thermal radiation model extends to the thermal design engi- 
neer the capability to define the inQuences of plume flow field and propellant 
composition on the magnitude and distribution of thermal energy emitted by 
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dual SRB plumes. Also, the capability of delineating the influence of surface 
shading on the magnitude of SRB radiant heating rates is provided. These 
capabilities are provided with a computer code that is easy to input for pur- 
poses of performing a detailed and extended parametric analysis of SRB 
thermal radiation in the base region of the Space Shuttle. 


1-4 


LOCKHEED ■ HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D496 763-1 


2. SRB PLUME STRUCTURE 

2.1 PARAMETERS INFLUENCING SRB FLOW FIELDS 

SRB exhaust plume flow fields are composed of the gaseous products of 
combustion as well as solid and/or liquid particles of Although the 

flow field of a single plume may be considered to be axisymmetric, calculation 
of the distribution of particle and gaseous properties in SRB plumes is a diffi- 
cult task because the properties of the gas and particles are not contiguous at 
any location. The particles are accelerated in the nozzle plume flow field due 
to pressure gradients and drag forces exerted on them by the gaseous flow and 
they retain temperatures that are generally higher than the gas temperature 
because of their thermal capacitance. Realistic flowfield calculations require 
accurate modeling of the interchange of momentum and thermal energy between 
the gaseous flow field and the dispersion of particles immersed in the 

gases in order to establish distributions of particles and gas properties in the 
plume. 

A number of parameters influence SRB plume flowfield structure. One 
of the more influential parameters is the ratio of pressure in the SRB com- 
bustion chamber to the ambient pressure (P^/P^). This ratio determines, 
to a large degreee, the amount of gaseous expansion that takes place in the 
plume and subsequently the plume's overall size. Analytical studies of SRB 
plume flow fields indicate that the initial expansion angle, y, (measured rela- 
tive to the plume axis) of the plume at the nozzle exit plane can be as large 
as 35 degrees. Figure 2-1 shows the effect of initial plume ex- 

nansion These data were obtained from several flowfield calculations using 

^ ( 1 ) 

Lockheed's two-phase plume flowfield computer code. At the lower values 

of P /P the SRB plume contains shock structures and associated localized 
c' oo’ ■ 1 

regions of high temperature gases and, to a lesser extent, Ai^O^ P3^^*ticles 

with locally higher temperature. These shocks tend to increase the magnitude 
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of tViermal radiation leaving the pluine. At high altitudes and corresponding 
high values of P /P , the plume expands considerably which reduces gas 
temperatures. Typical variations of chamber pressure during a nominal 
trajectory are shown in Table 2-1 as a function of time. Also shown are the 
ambient pressure and the ratio P^/P^, 

The influence that P^/P^ exerts on the plume flow field is often re- 
ferred to as "an altitude effect" since the magnitude of the pressure ratio 
is largely determined by changes in P^. However, there are other character- 
istics of SRB plumes that are due to variations in altitude as well as vehicle 
velocity. At low altitudes viscous mixing between the quiescent atmosphere 
and the plume leads to afterburning at the plume boundary. This occurs be- 
cause the plume tends to be fuel -rich due to incomplete burning of the pro- 
pellant in the combustion chamber. Afterburning caused by the mixing of 
oxygen in the atmosphere with the fuel-rich plume can greatly increase the 
heat load to the base plane of the Space Shuttle due to the plume boundaries 
proximity to the base region. Another parameter influencing plume structure 
that is somewhat dependent on altitude is the approach Mach number. The 
approaching freestream tends to lessen the amount of plume expansion at high 
values of freestream dynamic pressure. On typical trajectories freestream 
dynamic pressure is maximum at 35,000 to 45,000 feet and then diminishes 
with altitude. To include these altitude effects on the SRB flow field in the 
thermal radiation calculations, a sophisticated treatment of the flowfield 
problem is required. To a large extent, this is available in Lockheed's two- 
phase SRB flowfield code. Some SRB flowfield predictions from this code are 

shown in Section 2.3. 

2.2 GEOMETRY MODEL FOR SRB PLUME STRUCTURE 

The properties of the gaseous constituents and the Af^O^ particles vary 
along the axial and radial coordinates of the plume. To make these variations 
in plume properties easily handled by a math model for the solution of the 
radiation transport problem, it is necessary to represent the spatial varia- 
tions by means of relatively simple geometries but still provide a representa- 
tive description of the plume. This can be provided by considering the plume 
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I 
I 

I Table 2-1 

I TYPICAL VARIATIONS IN CHAMBER PRESSURE, AMBIENT PRESSURE 

AND P^/P_ DURING NOMINAL SPACE SHUTTLE TRAJECTORY 

V 00 


Time 

(sec) 

P 

00 

(psia) 

P 

o 

(psia) 

P /P 

2 

14.7 

768.8 

52 

4 

14.7 

781.5 

53 

6 

14.6 

792.1 

54 

10 

14.3 

802.7 

56 

16 

13.5 

814.1 

60 

18 

13.2 

817.6 

62 

20 

12.8 

812.2 

63 

22 

12.4 

774.5 

63 

24 

11.9 

744.9 

62 

32 

9.9 

661.1 

66 

40 

7.8 

593.0 

76 

48 

5.8 

543.7 

99 

50 

5.4 

534.3 

98 

52 

4.9 

525.0 

109 

54 

4.4 

527.9 

119 

56 

4.0 

532.0 

132 

58 

3.6 

534.5 

147 

64 

2.6 

544.7 

207 

72 

1.5 

568.5 

367 

80 

0.9 

584.3 

647 

88 

0.5 

596.3 

1190 

94 

0.3 

607.2 

2025 

96 

0.26 

611.2 

2350 

100 

0.18 

610.8 

3390 

1 10 

0.07 

570.0 

8140 

114 

0.05 

488.8 

9780 


2-4 


LOCKHEED ■ HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D496763-I 


to be constructed of homogeneous regions wherein appropriate values of gaseous 
and particle properties are defined. These homogeneous regions are stacked 
along the axis and radius of the axisymmetric flow field. A region is bounded 
along the radial coordinate by concentric, conical surfaces emanating from a 
common vertex and is bounded along the axial coordinate by planes at right 
angles to the conical surface axes (see Fig. 2-2). Properties are defined in a 


I 



Fig. 2-2 - Geometry Model for SRB Plume Structure 
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cylindrical coordinate system with its origin at the vertex of the cones. 
Average properties for each region are calculated based on assigned prop- 
erties on the four corners of each region. The assigned values for these 
properties may be determined from a two-phase flow calculation for the 
gasdynamic definition of the SRB plume using Lockheed -Huntsville's two- 
phase plume flowfield calculation. However, any valid definition of the 
spatial distribution of plume properties may be used as input data. 

The scheme used for defining properties for each homogeneous region 
of the plume makes use of volumetric averaging. One plume region is shown 
in Fig. 2-3. The property P(h, r) is taken to be a linear function of h and r 
as shown in Eq. (2.1). 



Fig. 2-3 - Property Definition for Plume Region 
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P(h, r) = P . + 

J Ah 


(h-H.) + ^ 
^ Ar 


(r - rp 


( 2 . 1 ) 


J 


where 


AP 

Ah 


P - P . 

1+ 1..1 i.j 

H. 1 - H. 
1+1 1 


( 2 . 2 ) 


and 


AP 

Ar 


fp. . , + AP 1 (h-H.)l - fp. . + AP . (h-H.)l 


Ar 


(2.3) 


The average value for the property P is found by integrating P(h, r) over the 
volume of the region and dividing by the volume of the region. 


P ^ y' P(h, r) dv/f dV (2.4) 

V V 

Average values for the following plume properties are assigned to each plume 
region by means of the above averaging technique; 

3 3 

• Particle number density, parts/ft (parts/cm ) 

• Particle temperature, R (K) 

• Particle radius, ft (microns) 

• Gas temperature, R (K) 

• Static pressure, atm 

• Mole fraction for CO, (-) 

• Mole fraction for (-) 

• Mole fraction for 1^2®’ 

• Mole fraction for HCi, (-) 
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The above properties result from a two-phase flowfield calculation for the 
SRB plume using Lockheed*s code. The mole fractions of the principle 
gaseous emitters (CO, ^^2’ ^2^ HCjp) are influenced by the method with 
which the plume chemistry is modeled. Available options for simulating 
plume chemistry are provided by using equilibrium chemistry properties or 
non -equilibrium chemical kinetics. However, representative calculations of 
SRB plume chemistry is provided by using equilibrium chemistry. 

2.3 EXAMPLES OF SRB FLOW FIELDS 

Some results of Lockheed*s two-phase flowfield calculations for SRB 
plumes are presented in this section. Data are shown for a sea level plume 
without an afterburning calculation, for a sea level plume with afterburning, 
and for a plume at an altitude of 72,000 feet without afterburning. Figures 
2-4, 2-5 and 2-6 show the boundaries for the gaseous and particle plumes. 

The particle plume boundary is defined by the limiting streamline of the 
particle trajectories at the edge of the plume. The limiting streamline (the 
streamline beyond which there are no more particles) is usually defined by 
the smaller particles in the flow field which are accelerated more rapidly in 
the radial direction by pressure gradients in the expanding plume than are the 
large particles. Defining the gaseous boundary is a much more difficult task, 
particularly for low altitude plumes where mixing with the ambient atmosphere 
occurs along with afterburning. The plume boundary for the sea level plume 
(Fig. 2-4) can be taken as the inviscid boundary when a mixing and afterburning 
calculation is not considered. The gaseous plume boundary for the afterburning 
simulation (Fig. 2-5) is taken to be the region where the gases in the afterburning 
shear layer drop below 3000 R (1667 K). The gaseous plume boundary for the 
plume at 72,000 feet (Fig. 2-6) is defined by the inviscid plume. 

Also shown in Figs. 2-4, 2-5 and 2-6 are the grids for approximating 
the property variations with homogeneous regions stacked along the axis and 
radius of the plume. Flowfield properties are supplied at the intersections 
of the conical rays and the horizontal lines. For grid points that lie outside 
the plume boundaries, the property values are simply set equal to zero. 
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Fig. 2-5 - Gaseous and Particle Plume Boundaries for SRB 
Sea Level Plume with Afterburning 
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This is particularly evident for the 7Z,000 foot plume where the large initial 
expansion angle at the nozzle exit plane requires a large cone angle for the 
outer grid rays. However, the plume flow field expansion is retarded down- 
stream (at Z/R^^ = 6) by the approaching freestream which has a Mach num- 
ber of 2.5 at the 72,000 foot trajectory point. 

Figure 2-7 shows the radial distributions of the properties of the 
particles that enter into thermal radiation calculations, namely particle number 

density (N), particle radius (r ) and particle temperature (T ). These data 

P p 

are shown at three axial stations in the plume Z/R^^ = 0.8, 4.0 and 11.0. The 
particle radius tends to drop at the edge of the plume because smaller prop- 
erties "migrate” to the edge more rapidly than do larger particles as the plume 
expands down its axis. The particle temperature drops at the edge due to the 
presence of smaller particles which more quickly adopt the cooler gas tem- 
perature. The particle number density diminishes as the flow proceeds down 
the plume axis because of the expanding volume of the plume. The particle 
distributions shown in Fig. 2-7 are used for both sea level plumes calculated 
for afterburning and without afterburning. 

The gas temperature and pressure and the mole fractions for the four 
principle gaseous emitters (CO, ^^ 2 ' ^2^ and HCf) are shown in Figs. 2-8 
and 2-9 for the sea level plume without an afterburning calculation and with 
an afterburning calculation. The gas temperature and the mole fractions are 
strongly influenced by the afterburning calculation as is evident by comparing 
Figs. 2-8 and 2-9- If mixing in the shear layer between the ambient atmos- 
phere and the plume flow is not considered, the combustibles in the plume are 
maintained in a chemically frozen state as is indicated by the lack of change 
(with respect to r/R^^) in the mole fractions shown in Fig. 2-8. However, if 
afterburning is simulated, a strong combustion process occurs at the plume 
boundary as is evident from the large changes in the mole fractions at the 
edge of the plume that results as CO and HC^ burn to CO^ and H^O (see Fig. 
2-9). The gas temperature is raised at the edge of the plume. Also, the 
volume of hot gases that will radiate back to the base plane is increased 
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Fig. 2-7 - Radial Particle Property Distributions for Sea Level Plume 
Without Afterburning 
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because the atmosphere at the edge of the plume enters into the combustion 
process. Comparison of Figs. 2-8 and 2-9 shows that the afterburning calcu- 
lation causes the effective plume to spread considerably at axial stations re- 
moved from the nozzle exit plane (Z/R^^ = 11,2), The gaseous plume boundary 
is taken to be at the region where the gas temperature drops below 1667 K, 

Property data for a SRB plume at 72,000 feet are shown in Fig. 2-10 
{for particles) and Fig, 2-11 (for gases). The 72,000 feet plume has no after- 
burning simulation as the ambient atmosphere is rarefied to the point where the 
amount of afterburning is expected to be low. 
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3. SPACE SHUTTLE BASE HEATING FROM DUAL SRB PLUMES 

The two SRBs astride the external taink (ET) on the Space Shuttle con- 
tribute significantly to the heat loads incident on spacecraft surfaces in the 
base region. In order to calculate radiant heat rates to these surfaces with 
realistic representation of the SRB plume radiation field, it is necessary as 
a first priority to establish the magnitudes and distribution of flowfield prop- 
erties which contribute to the radiant energy leaving the plume. This infor- 
mation can be provided by Lockheed's two-phase flowfield code as described 
in Section 2. 

3.1 DUAL SRB PLUME GEOMETRY MODEL FOR THERMAL RADIATION 

Two SRB plumes are considered in the SRB thermal radiation model. 

The plume locations are defined in a central Cartesian coordinate system 
centered between the dual SRB nozzles. The X^ coordinate is 
parallel to the plume Z-axis at zero gimbal. The X^ coordinate passes 
through the center of the exit plane of the two SRB nozzles and the coordi- 
nate is at right angles to X^ and X^. Both plumes are considered to be axi- 
symmetric and the distributions of plume properties along the axial and radial 
coordinates are approximated by defining homogeneous regions along the axis 
and radius of the plume as outlined in Section 2. The geometry model allows 
the SRB plumes to be gimbaled parallel to one another in the yaw and pitch 
planes or they may be gimbaled independently with different yaw and pitch 
gimbal angles. Figure 3-1 shows the dual plume geometry model for zero 
gimbal, parallel gimbal and independent gimbal in the yaw plane. The same 
type of gimbal angle options are available for gimbal in the pitch plane. Simul- 
taneous gimbal in both yaw and pitch planes may also be modeled. 

The SRB plume may be considered to extend several hundred feet down- 
stream of the nozzle exit plane. However, for purposes of calculating radiant 
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heat loads to the base region of the Space Shuttle, only a finite portion of the 
plume need be considered. Often the nature of the SRB flow field can be used 
to determine how much of the plume should be used in calculating radiant heat 
transfer rates. The variation of plume temperature along the axial coordinate 
provides a good indication of how significantly portions of the plume will con- 
tribute to thermal radiation in the base plane. Usually a plume axial length of 
10 to 15 exit radii is sufficiently long to provide good estimates of radiant heat- 
ing. Longer plume lengths are not necessary because of the diminishing view 
of the plane to the base region. As the plume length considered for heating 
rate calculations increases, it is likely that the dual plumes used in the heat 
rate calculations will intersect. Plume intersection will certainly occur for 
the finite plume lengths used to simulate plume structures of high altitude 
plumes because of the high degree of plume expansion at altitude. The dual 
SRB plume geometry model considers the possibility of plume intersections 
as will be discussed later in this report. 

3.2 THERMAL RADIATION FOR GRAY MEDIA 

The spatial distribution of radiant energy leaving the SRB plumes is 
established by the Monte Carlo method. This method considers the radiant 
energy to be divided into several discrete bundles of equal energy. These 
bundles are considered to obey the physical laws that govern the emission, 
absorption, and scattering of photons. In the Monte Carlo method a large 
number of these energy bundles (referred to as the sample size) are traced 
as they undergo simulated emission, scattering, or absorption and subsequent 
re-emission until they eventually escape from the boundaries of the plume. 

Upon escaping the plume, the energy bundle travels along a straight line (re- 
ferred to as a line of sight) as would a photon leaving the plume. This line 
of sight is easily defined by a linear algebraic expression. In the computer 
simulation of tracing the energy bundles as they leave the plume, the alge- 
braic definition of the line of sight may be checked to determine if the energy 
bundle intersects one of the simulated surface geometries that represent 
spacecraft components. If the bundle strikes one of the Space Shuttle surfaces, 
the energy contained in the bundle is considered to contribute to the radiant heat 
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load on that surface. The problem of spacecraft surface heating will be 
discussed later. This section addresses the formulations required to 
model magnitude and distribution of thermal radiation leaving SRB plumes. 


3.2.1 Thermal Radiation from a Participating Medium 


An expression for the magnitude of radiant energy emitted locally in 
SRB plumes may be obtained by considering the attenuation of radiation 
intensity in the plume. Consider a volume element of the SRB plume. Con- 
tained within the volume element are a homogeneous dispersion of Af O 

2 3 

particles immersed in the gaseous products of combustion. Although there 
are several gaseous constituents present, only CO, CO^, H^O and HCf con- 
tribute significantly to thermal radiation. Radiation intensity at a given 
wavelength in the volume element will be attenuated in a manner represented 
by the following expression 


5 . 

dx 


- •'a, X + X> h 


(3.1) 


where 

N = particle number density 

^ = spectral absorption cross section of the particles, and 
^ = volumetric absorption coefficient of the gas 


Imposing the constraints of local thermodynamic equilibrium (LTE) requires 

that the net change of radiation intensity be zero along a given line of sight. 

This condition is satisfied when there is an amount of energy emitted by the 

volume element equal to that absorbed. Let J represent local emission in 

A 

the volume element into the solid angle about the line of sight. Then imposing 
LTE requires 


^ "a, X + 'a ° 


(3.2) 


or 


"x = "a, X + ^a, x> 'x 


(3.3) 
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Equation (3.3) is applicable for an incremental solid angle, dco , about the 
line of sight under consideration.. Integrating Eq. (3.3) over 4n steradians 
and over a finite volume gives 


J = 4 (N CT E^ n -X ^ 

X ' a, X b, X a, X b, X 


(3.4) 


where 


and E® = spectral black body functions for the 
^ ^ particles and gases, and 

V = volume under consideration. 


Equation (3.4) represents the spectral radiant energy emitted by a 
homogeneous volume of particles and gases under conditions of LTE. 

Rigorous calculations of thermal radiation leaving a volume of emitting, 
absorbing, and scattering medium would account for the influence of wave- 
length dependent properties. However, for purposes of calculating the total 
thermal radiation heat load, the gray approximation to the radiation trans- 
port phenomenon is often sufficient. Simulating the coupled radiation from 
the particles and the gaseous emitters in SRB plumes as a gray medium re- 
quires that appropriate wavelength averaged values be obtained for the ab- 
sorbing and scattering properties of the medium. The gaseous constituents 
of the plume are almost pure absorbers of radiant energy. Conversely, the 
Ai particles are almost pure scatterers of radiant energy when they are 
in the solid phase. Liquid ^{^^3 however, can have significantly 

large absorption cross sections. 

3.2,2 Gray Absorption Coefficients for SRB Plume Gases 

Gases tend to absorb and emit radiant energy in discrete wavelength 
bands. The absorbing properties of the gases can change by several orders 
of magnitude over a small wavelength interval. An appropriate averaging of 
absorbing properties over the wavelength spectrum must be performed 
order to calculate meaningful radiant heating rates. Abu-Romia and Tien 
have made a study of appropriate mean absorption coefficients to be used in 
thermal radiation calculations for gases. It is recommended in Ref. 2 that the 
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Planck mean absorption coefficient be used for gases that approach the opti- 
cally thin limit. Considering the gaseous phases alone, the thin limit is 
approximated in SRB plumes and the Planck mean absorption coefficient be- 
comes a valid means of removing spectral dependence from plume properties 
in the thermal radiation calculations. 

The Planck mean absorption coefficient is defined as 

00 

/ ^a,X ^b, A T) 

K^(T,P) =-2 35 (3.5) 

/ E,,^a.T)dx 

o 

wherein the spectral variations of the volumetric absorption coefficient is 
weighted by the blackbody function. The Planck mean absorption coefficient 
will be independent of wavelength but will still be a function of the partial 
pressure and temperature of the gas under consideration. Figure 3-2, taken 
from Ref. 2, presents the Planck mean absorption coefficients for CO, CO^ 
and H^O as functions of gas temperature for a partial pressure of one atmos- 
phere. The Planck mean absorption coefficients shown in Fig. 3-2 are used to 
represent CO, CO^ and H^O as gray gases in the thermal radiation model. 

Similar data for HCi are needed but are not available. In order to include 
HCI in the thermal radiation calculations, the Planck mean absorption coef- 
ficient for HCf is taken to be that of CO. This approximation should provide 
a preliminary estimate of the absorption coefficient for CO in that the spectral 
absorption coefficient for HCf is lower than that for CO but, on the other hand, 
the band for HCi subtends a greater fraction of the blackbody function than does 
CO. These two effects should mitigate against each other and make the approxi 
mation more valid. 

The Planck mean absorption coefficients shown in Fig. 3-2 were 
calculated for a partial pressure of 1 atm. In order to use these coefficients 
directly in calculations of thermal radiation, it is necessary to ratio the coef- 
ficient in Fig. 3-2 from 1 atm to the partial pressure of the gas being modeled. 
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Taking this into consideration, the absorption coefficients for CO, CO^, H^O 
and HCf are added linearly in order to obtain the absorption coefficient for the 
entire mixture as in Eq, (3.6). 


K 

mixture 


P 

1 atm 



M.F.. 


1 


(3.6) 


where 

P = static pressure of the gas mixture 

K = linear absorption coefficient of constituent gases 
^i HCj? , and 

M.F. ^ = mole fraction of constituent gases. 


3.2.3 Scattering and Absorbing Properties of Particles 

Af^O^ particles can possess substantial absorption characteristics as 
well, depending on the magnitude of the imaginary portion of the index of re- 
fraction. For liquid particles of the index of refraction can be large 

which leads to relatively large absorption cross sections. On the other hand, 
solid At^O^ particles are almost pure scatterers of radiant energy. The cal- 
culation of the magnitude of the absorption and scattering of electromagnetic 
radiation by particulate matter, such as the distribution of Aft^O^ particles in 
the SRB plume, can be done by using the Mie theory; its application to calcu- 
lating the absorbing and scattering properties of Al^O^ clouds is contained in 
the following paragraphs. 


The problem in scattering of electromagnetic radiation by particulate 
matter, such as spherical particles immersed in a vacuum, is to relate the 
properties of the particles (i.e., size, shape, index of refraction) to the spatial 
distribution of the scattered radiation and to the absorption of radiation by the 
particulate matter. Basically, there are three regimes for the solution of this 
problem. The limits of each regime are dependent on a size parameter, namely, 
the ratio of the circumference of the scattering particle to the wave length of 
the incident radiation 

2it r 


K 


a 



(3.7) 
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For wave lengths of incident radiation orders of magnitude larger than the 
dimensions of the particle, the Rayleigh approximation to scattering of the 
radiation is applicable. For wave lengths orders of magnitude less than the 
dimensions of the particle, classical optics may be applied to the scattering 
of the incident electromagnetic radiation. When the wave length of the incident 
radiation and dimensions of the scattering particle are on the same order of 
magnitude, recourse to sophisticated treatment of the scattering and absorption 
phenomenon at the particle is required. 

The Mie theory treats the case where the wave length of incident radia- 
tion and the size of the particle are comparable. Mie theory is concerned 
with a very idealized situation, that is. a sphere composed of homogeneous, 
isotropic material embedded in a homogeneous, isotropic, dielectric, infinite 
medium irradiated by a plane wave propagating in a specified direction. The 
dielectric properties- of the matrix medium imply that this medium is non- 
absorbing. The gaseous constituents in the SRB plume are physically repre- 
sentative of such a medium. For a sphere composed of a conducting material, 
that is, a non -dielectric material, any incident radiation on the sphere is partly 
absorbed, scattered and transmitted by the sphere. The Mie solution predicts 
the scattering and absorption cross sections (or efficiencies) as well as the 
spatial distribution of the scattered intensity. The efficiencies are the ex- 
tinction, absorption, and scattering cross sections normalized by the geom- 
etric cross section 

Q .3 ^e, a, scat (3 8 ) 

''^e, a, scat ^ ^2 

P 

The Mie solution is effected in the form of complicated series that involve the 
Riccati-Bessel functions and Riccati-Hankel functions of increasing order. 

A computer program for the Mie solution is discussed and listed in Ref. 3, 

The extinction and absorption efficiencies and the absorption to extinction 
ratio for particles are shown in Fig. 3-3 as a function of wavelength. 

Particle radii of 5, 6 and 7 microns are considered and the index of refraction 
for is taken to be 1.8 + .008i, The magnitude of and as predicted 

by Mie theory are actually a function of one variable, namely - Zir r^/X, but 

3-9 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 





i 

s 


Fig. 3-3 - Spectral Variation of Extinction Efficiency and the Absorption/Extinction Ratio 
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the physical influence of and X are perhaps best pointed out in parametric 
form as in Fig. 3-3. It is clear from Fig. 3-3 that the (and therefore a^) 
is influenced substantially by wavelength of the radiation field. Although the 
wavelength dependence of cr^ and for At. 2 ^ 0 ^ is less critical than the de- 
pendence of linear absorption coefficients for gases, a mean value for and 
o for the particle dispersion should be used. The Planck mean is again used 

for this purpose as shown in Eq. (3.9); 

00 




aa(X),a^(X) E^(X. T) dX 
00 

J E^iX, T) dX 


(3.9) 


O 


3.2.4 Total Thermal Radiation Emitted by the SRB Plume 


Using the Planck mean values for the absorbing properties of the con- 
stituent gases of the SRB plume and the At^O^ particles, an estimate may be 
made of the spontaneous emission of radiant energy of the participating medium 
that comprises the SRB plume. Using Eqs. (3.6) and (3.9) an expression may be 
obtained for the gray approximation for radiant energy leaving one region of the 
plume as shown in Eq.(3.10): 

J = 4 (N a E^ + K E®) V (3.10) 

^ a b a b 


The total radiation emitted by the entire volume of the plume is taken to be the 
sum of the radiant energy emitted by all regions of the plume. 

all 

regions 

^total = “ Z) (n % 
i = l 


+ K 


^b 


V. 


(3.11) 


3.3 APPLICATION OF THE MONTE CARLO METHOD TO SRB PLUME 
THERMAL RADIATION 

The Monte Carlo method is generally used for solving problems in 
which the physics of the phenomenon of interest is understood but in which 
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the equations that describe the physics are difficult to solve due to complicated 
boundary conditions* Application of the Monte Carlo method to thermal radia- 
tive transfer problems involves establishing a math model of the physics that 
govern the thermal radiation processes* The physics modeled are those that 
characterize the emission, scattering and absorption of photons. Basically, 
the solution consists of a group of a large number of photons throughout the 
system. For an emitting, scattering, and absorbing medium such as a solid 
rocket plume, the groups of photons, herein referred to as energy bundles, 
are considered to be emitted at appropriate points throughout the plume. The 
bundles are then traced along probable paths that may lead to several scatter- 
ing and/or absorption events with the particles and the gaseous con- 

stituents of the plume. An absorbing collision will, in effect, terminate 
interest in the subject bundle in that the bundle cannot escape the plume and 
there is no possibility that the energy in the bundle will be imparted to a target 
of interest. When an energy bundle is absorbed, a second bundle is considered 
to be emitted at the absorption site in order to satisfy LTE, If an energy bundle 
escapes the plume either directly upon emission or after several scatterings 
and r e -emis sions , its path upon escaping the plume is analyzed to determine 
if it hits or misses some predetermined target. If the bundle strikes some 
target which models a Space Shuttle surface of interest, a portion or all of the 
energy in the bundle may be considered to be absorbed by the target surface 
depending on the manner in which the absorbing properties of the surface are 
simulated. 

In order to perform a Monte Carlo simulation of the thermal radiation 
leaving SRB plumes it is necessary to properly select emission sites and 
directions as well as to set up geometry models to track energy bundles 
within the plumes. This section addresses these topics and presents ex- 
pressions for calculating heat flux to target surfaces. 

3.3.1 Distribution Functions 

The emission, scattering and absorption (and re-emission) of energy 
bundles must be simulated in accordance with the physical laws known to 
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govern photons. Howell demonstrated that appropriate selection of random 
events may be achieved by setting the cumulative distribution function for the 
variable equal to a set of random numbers that is uniform between 0 and 1 
and then solving for the variable. Sets of random numbers are easily gen- 
erated by a digital computer. Random number generators are usually sup- 
plied in subroutine form as part of the digital computer software. The 
nature of the random number generator available on the MSFC Univac 1108 
is reviewed here. A more detailed discussion is available in Ref. 5. The 
random number is generated by squaring a six-digit integer and then taking 
the middle six digits from the squared value which is then divided by 10 re- 
sulting in a number between 0 and 1. The six digit number which initializes 
the random number array is supplied as an input value (NSTART) for the 
Univac random number generator. The same sequence of random numbers 
is generated if the same value for NSTART is used as an initializing value. 

Such a sequence is referred to as pseudo -random and the process of the multi- 
plicative congrucntial method. The method supplies a distribution of random 
numbers which is uniform over the interval from 0 to 1 . 

So, a uniform set of random numbers is conveniently supplied as part 
of the MSFC Univac 1108 software. Now, if the cumulative distributions for 
the variables of interest in solving the radiation transport problem can be 
defined, the radiation transport problem can be solved with the Monte Carlo 
technique. References 4 and 6 outline the method for obtaining cumulative 
distribution functions which amounts to integrating the probability density 
function as 

X 

F(x) = / P(x) dx (3.12) 

-00 

where 

P(x) = probability density function 

F(x) = cumulative distribution function, and 

X = variable of integration. 

P(x) defines the probability that an event will lie between x and x + dx. As 
such, the integral of P(x) over all x is unity. Thus, the cumulative distribution 
function has a value between 0 and 1 for any value of x. 
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Cumulative distribution functions must be derived which simulate the 
following phenomena in order to calculate the distribution thermal radiation 
leaving SRB plumes: 

• Uniform distribution of emission sites over each homogeneous 
region of the plume 

• Isotropic emission from a volume element 

• Isotropic scattering 

• Anisotropic scattering, and 

• Attenuation of radiant intensity in a homogeneous medium. 


Selection of Ernission Sites 


Each region of the plume emits an amount of energy as thermal radia- 
tion as defined by Eq, (3,10). The number of spontaneous emissions that are 
allowed to take place in each region is determined by the fraction of the total 
energy leaving the plume that is emitted by the given region as shown in Eq. 
(3.13): 



J. 

1 


I N CT E?i 


+ K E' 


IV. 


e, 1 


all 

r eg i on s 

E 

i = l 


(3.13) 


N o- Ef; 
{ a b 


+ K E^ 
a b 


V. 


where ^ is the number of energy bundles emitted by region i. In this 
manner the number of energy bundles emitted by each region is made pro- 
portional to the amount of radiant energy leaving the region. 


It is also necessary that a uniform distribution of emission sites is 
achieved within a given region. The SRB plumes are assumed to be axi- 
symmetric and their structure is defined in a cylindrical coordinate system 
centered at the vertex of the conical surface boundaries. So, proper values 
must be assigned to the cylindrical coordinates r, ^ and h, A uniform dis- 
tribution within the homogeneous region implies that 
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^ dV r dr dh 

N " V ~ V 

e 


(3.14) 


The volume of a region of the plume may be expressed as 
Zt -H. , , ^h tany. 

C C I ’ 

- i } ' 

o H. 


j+1 2 2 + 

r dr dh = jr (tan Yj) 3 (3.15) 


h tany^ 


Reference to Fig. 3-4 will aid in defining the variables. Now 
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2r dr 


h (tan y 


j+1 


- tan 




7f 


3h^ dh 


4l 


1 


(3.16) 


The fact that the coordinates r, 4> and h are independent and that Eq. (3.16) 
integrates to unity, implies that Eq. (3.16) is the joint probability density 
function for the selection of emission sites within the prescribed volume. 

The probability functions for each coordinate may be taken as (since coordi- 
nates are independent) 

F(*) = / = 2.U^ (3.17) 


F(h) = 


n 
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H. 


, = u. =» h 
1 + 1 1 




^ - H?) 
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+ 1 ■ ""i ' "h 


(3.18) 


F(r) = 
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/ 


2r dr 
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h^ (tan^y^^j - tan*'yj) 
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(3.19) 
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Equations (3,17), (3.18) and (3.19) have been attained by integrating (with re- 
spect to the coordinate) the probability density function for each coordinate 
in order to obtain the cumulative distribution function. The cumulative dis- 
tribution function is then set equal to a random number (U^, U^, U^). Using 
these equations, the random selection of emission sites may be put into effect 
with assurance that the emission sites so selected will be uniformly distributed 
over the given volume, 

• Isotropic Emission and Scattering 

Once an emission site has been selected it is necessary to select the 
direction of travel of the subject energy bundle defined in a spherical co- 
ordinate system centered at the event site. Isotropic implies that the direction 
selected will have no bias for any given direction. This requires that the prob- 
ability of emission into any incremental solid angle, du) , be identical. Thus, 



Equations (3.21) and (3.22) govern the selection of direction coordinates for 
isotropic emission as well as isotropic scattering since the selection of 
coordinates for isotropic scattering is governed by the lack of direction bias. 


• Anisotropic Scattering 

A Mie theory prediction of the scatter function S(r], 9) for anisotropic 
scattering from particles is shown in Fig. 3-5. A scatter function for 
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isotropic scattering is also shown for comparison. The scatter function is 
defined as the probability that the scatter direction relative to the direction 
of the incident radiation will fall into the solid angle, du)'. Thus, 

S(n>9)f7 = ^ S{rj,6) sinrj' drj' d0- 

If unpolarized radiation is considered, an even distribution of scattered radia- 
tion about the incident direction results. Then, 

h' 

F(r)') = S(rj-) sinr,' dp' = U^. (3.23) 

o. 

and 

F(0') = 2^ y de' = Uj, =:* 9’ = (3.24) 

o 

The anisotropic. nature of scattering by Al^O^ is manifested principally 
by a strong bias for scattering into the solid angle in the forward direction. 

In order to simulate this phenomenon some means of solving the integral ex- 
pression in Eq. (3.23) must be performed. Stockham^^^ has discussed a con- 
venient method for doing this. The area under the curve of S(p') sinp' versus 
p' shown in Fig. 3-5 is normalized to unity. This area may be represented by 
the sum of the areas of triangles and rectangles which constitute the total area 
under the curve. Then, the polar angle for anisotropic scattering may be ob- 
tained from the expression 

A(p') = U^, 

where 

A(p') is the functional relationship for the summation of 
triangular and rectangular areas under the curve of 
S(p') sinp' versus p'. 

• Attenuation of Radiation Intensity 

The intensity of a beam of radiation traversing an absorbing and scatter 
ing medium is attenuated in proportion to the intensity of the beam and the 
optical properties of the medium. This is expressed by Eq. (3.1) wherein 
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spectral values are considered. Rewriting Eq. (3.1) for a grey medium com- 
posed of AH^O^ particles and portions of CO, CO^, and HCi gives 


dS 


(N + K^) I 


or 


(3.25) 


I 


(No + K ) dS 

St 3. 


(3.26) 


In a homogeneous region, N, a , and K are constant and 


/ f = - (N + K^) y dS 

r o 


(3,27) 


which integrates to 


I ■ ^ 

o 


(3.28) 


Equation (3.28) represents the fraction of the original beam intensity 
at Station S relative to the beam intensity at Station 0. The fraction of 
intensity removed prior to S is 

'-f = 

o 

and Eq. (3.29) can be thought of as the probability that the path length will be 
less than or equal to S and thus constitutes the cumulative distribution func- 
tion for defining S. 

- (N a + K ) S , 

F(S) = l- e ^ ^ =U„=^S'= -^ 1 = 3 !n U„ (3.30) 

^ P N CT + K ^ 

a a 

Since Ug is equivalent to (Ug - 1) in the random sense, S^ represents the 
probable length of travel in the medium. 


1 - e 


-(No + K ) S 
a a 


(3.29) 
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It is convenient to normalize the length dimensions with the exit radius 
of the SRB nozzle in order to more conveniently represent the optical prop- 
erties of the plume in terms of optical depths (really optical radius) which are 
dimensionless. Dividing Eq. (3.30) by 

S = !n Ug (3.31) 

P «ex (N + K^) 

The optical properties of the A£^0^ particles and gaseous constituents can 
now be expressed in terms of optical depths based on the nozzle exit radius, 

R . 

ex 


1 

' ■ 


T 

= N (7 R 

(3.32) 


P 

a ex 


I 


T 

= K R 

(3.33) 



g 

a ex 


1 

a 


T = Tp 4- Tg 

= (N a + K ) R 
'a a ex 

(3.34) 


and 




: 1 


s 

= - - in Uc 

(3.35) 

- 


P 

T S 



• I Equation (3.35) defines, in terms of random numbers, the probable 

[ length of travel of an energy bundle prior to either a scatter or absorption 

^ I event in a single homogeneous region. However, the SRB plume is not ho- 

I mogeneous and is modeled as a composite of homogeneous regions as dis- 

■ cussed in Section 2. The probable length of travel for an energy bundle may 

take it across several different homogeneous regions before it either suffers 
• a scatter or absorption event or escapes the plume. Each region of the plume 

® has its own unique optical properties and attenuates the radiation intensity 

. passing through it at different rates. Consider a beam of radiation passing 

I through a series of homogeneous media as sketched in Fig. 3-6. The line of 

sight passes through five different regions numbered successively 1 through 5. 
I The length of travel in each region and the optical properties of the region 

define the rate of attenuation of radiant energy from the beam. 
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Fig. 3-6 - Line of Sight of Energy Bundle Traversing Several Plume Regions 


The attenuation from S = 0 to any station along the line of sight is 
expressed as 

^1 '’'2 ^2 ’’^i ^i^ 

I " ^ 
o 


(3.36) 


The probability of collision along the line of sight is 


1 - — = 1 - e 


-(Ti +T 2 +[Sp-(Si +S2+...+S^_j)] T^) 


(3.37) 


= F(S) = U, 
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or 


- Sn (Ug) 



(3.38) 


where the subscript n refers to the region in which the scatter or absorption 
event occurs. Solving explicitly for the probable length of travel 


n-1 

S = V s 

p jLj 1 

i = l 


n-1 


i = l 


n 


(3.39) 


It is again noted that all lengths are normalized with respect to the exit 
radius . 


• Selection of Scatter or Absorption Event 


As the life history of an energy bundle is followed during its residence 
within the boundaries of the plume, the bundle may be subjected to several 
scatter and absorption events. Once it is determined that the subject energy 
bundle will suffer a collision event, it is necessary to specify whether the 
event is an absorption or a scatter event. This is done by establishing the 
absorption to extinction ratio for the various regions of the plume. The Ai 203 
particles have a significant scattering cross section while the gaseous con- 
stituents of the SRB plume are pure absorbers. So an expression for the ab- 
sorption to extinction ratio for a homogeneous region may be expressed as 


A/E 


G /a + T„ 
e 'p K 


(3.40) 


A/E is compared with a random number, U. Tf U is less than A/E, the event 
is considered to be an absorption event. For U larger than A/E, the event 
is a scatter event. 
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3.4 TRACING ENERGY BUNDLES WITHIN THE PLUME 


It is necessary to track the location of energy bundles as they are 
emitted, scattered or absorbed and reemitted during their lifetime within 
the plume. It is also necessary to identify the region of the plume that the 
bundle is traversing in order to attenuate the radiation intensity at the appro- 
priate rate. The first step in this procedure is to calculate the elevations 
of intersections of the line of sight with the conical surfaces which define 
the boundaries between regions along the radial coordinate of the plxime. 


3.4,1 Calculation of Line-of-Sight Intersections with Conical Region 
Boundaries 

The method used to calculate the elevations along the Z-axis of inter- 
sections of the line of sight with the conical boundaries follows that used by 
Tien and Abu~Romia: ^ As a first step the distance is calculated between the 
axis of the conical surfaces and the point where the line of sight pierces a 
plane at right angles to the axis and passes through the cone vertex. Figure 
3-7 shows this intersection point, P^; li and locate the point within the 
plume where an emission, scatter or absorption/reemission event occurs. 
Using the law of cosines, the distance R between the plume axis and the 
piercing point is 

R = ^r^ + (h tanr^)^ ” 1 ^ cos(jr - (9 ~ ^ (3.41) 


Now referring to Fig. 3-8 the following geometric identities are 
apparent 


tany 



(3.42) 


tanr] 



(3,43) 


tanr^’ 



R +c 



(3,44) 
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Fig. 3-7 - Orthographic Projections of Line of Sight to Piercing Point 
in Plane Through Vertex 
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I 


2 



Fig. 3-8 - The Path Length Through a Conical Body 


3-26 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D496763-I 


Tlirouj^h use of Pythagoras theorem and the quadratic equation, the following 
oxpressions are obtained which give the Z elevations for intersections of the 
line of sight with the conical boundaries relative to a plane through the vertex 
in terms of y, rj, tq' and R. 


Z 


il 


Z 


i2 


Vtanv + R tann* 

r 2 2 2 

- (tany + R tann') + (1 - R ) (tan q 

- tan^y)] 


2 , 2 



tan q - tan y 


‘\/tany + R tan rt ' 

2 2 2 

+ (tanv + R tann’) + (1 -R ) (tan ri 

- tan^y)] 

_ 

2 2 



. 2 ^ 2 
tan t] - tan y 


(3.45) 


(3.46) 


Equations (3.4 5) and (3.46) apply only when q > Y • If rj < Y , the line of 
sight exists through the base of the plume. For r) < y the Z elevation of the 
intersection of the line of sight with the cone region boundary is simply the 
elevation of the base. If q > 90 deg, Eq. (3.45) is used to calculate the Z co- 
ordinate of the intersection. If y < q < 90 deg, Eq. (3.46) is used. These 
equations are used to calculate the intersection of the line of sight with all 
conical region boundaries used to construct the plume. 

3.4,2 Attenuation Along the Line of Sight 

Once an energy bundle is emitted, scattered or absorbed/reemitted, 
the line of sight that the bundle follows may cross several regions before 
leaving the plume. In order to calculate the probable length of travel as per 
Eq.(3.39), the length of travel (S^) in each plume region must be calculated. 
Figure 3-9 will aid in visualizing the manner in which is calculated. In 
the following discussion the planes at right angles to the plume axis used to 
bound plume regions along the axis are referred to as I surfaces. The 
conical boundaries are called J surfaces. An event site is shown in Region 
j - 2, J = 2. For ease in visualizing the method, the trajectories of the 
energy bundle are shown for the case where the line of sight passes through 
the plume axis (which corresponds to an azimuth angle of 9 = 0 or 180 deg). 
This expedient makes the lines of sight shown in Fig. 3-9 appear in true 
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Fig. 3-9 - Line-of-Sight Intersections with Plume Region Boundaries 
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length. However, the equations to be derived for calculating the length of 
travel in each region are valid for any 9 and ^ . 


The location of event sites in the plume are tracked in a cylindrical 
coordinate system centered at the common vertex of the conical surfaces 
used to construct the plume. The true length of travel in any region may be 
determined by identifying the Z -coordinate of the inter sections of the line of 
sight with the boundaries of the region and dividing the difference in elevation 
by the cosine of the polar angle, r], of the line of sight. These intersections 
may occur on any combination of the I and J surfaces. The possible combi- 
nations of intersections are with I-I surfaces, I- J surfaces, J-I surfaces or 
J-J surfaces. These combinations of line -of-sight surface intersections are 
shown in Fig. 3-9. However, note that there are only two combinations of 
line-surface intersections from the event site to the region boundary, namely 
site to I-surface or site to J-surface. The Z coordinate of the intersections 
with the I surfaces are available from data used to construct the plume struc- 
ture, namely in the array of HZ^ values. The intersections with the J sur- 
faces are available from either Eq. (3.45) or (3,46), The true length of travel 
in any plume region may now be conveniently calculated from 


S. 


Z ^1 - Z 
n-H n 

cos n 


(3.47) 


These values of may be used in conjunction with Eq. (3.39) in order to cal- 
culate the probable length of travel, S , of an energy bundle. 


3.4.3 Location of Scatter and Absorption Event Sites 

In the simulation of SRB thermal radiation using the Monte Carlo method, 
the energy bundles are tracked in a cylindrical coordinate system centered at 
the vertex of the conical boundary surfaces with the positive Z-coordinate 
aligned with the cone axis. The direction of travel of an energy bundle is 
assigned in a spherical coordinate system centered at the event site. Emis- 
sion sites within the plume are selected based on Eqs. (3.13), (3.17), (3.18) 
and (3.19). 
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Equation (3.13) is used to select the plume region in which the emission 
site occurs. Equation (3.17) selects the azimuth angle measured counter- 
clockwise with respect to a line parallel to the positive coordinate of the 
principal coordinate system and at right angles to the plume axis. The evalu- 
ation of the emission site, h, relative to the vertex of the conical boundaries 
along the Z-axis of a single SRB plume is determined by Eq. (3.18). The 
radial distance from the plume axis, r, is determined by Eq.(3.19). Once 
an emission site is selected, the energy bundle is followed along a line of 
sight with polar and azimuth angles defined by Eqs. (3.21) and (3.22). The 
polar angle, rj, is measured with respect to the positive Z-axis of the plume. 
The azim.uth angle, 0, of the line of sight is measured with respect to a line 
parallel to X ^ -coordinate and at right angles to the plume axis. A probable 
length of travel is then selected based on Eq. (3,39). The energy bundle is 
tracked along the line of sight until it either escapes the plume or until a 
collision event occurs. The occurrence of a collision event is determined 
when the predetermined value for is less than for any region. This 
circumstance implies that the energy bundle failed to escape the plume along 
the specified line of sight. Once a collision event occurs it is necessary to 
calculate the cylindrical coordinates of the new event site. Reference to 
Fig. 3-10 should be made for the following discussion. Figure 3- 10 is a 
schematic of the geometry involved in tracing energy bundles. The plume 
shown in Fig. 3- 10 is composed of a single region. However, the relations 
derived for tracking energy bundles from site to site are valid for a plume 
composed of any number of regions. 

The Z-coordinate of the new event site is obtained by adding the pro- 
jection of Sp onto, the Z-axis of the plume 


^2 “ h j + S cos r? 


(3.48) 


The new value for the r -coordinate is obtained from the law of cosines. 

Consider the r^, the new value of {r^ in Fig. 3-10), and the projection 

of S onto a plane at right angles to the Z-axis. The projection of S onto 
P P 
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such a plane is simply S sin 7 ^ and the law of cosines gives 

P 


= + (Sp sinr])^ - 2r ^ sinr; cos( 7 r- (9 


(3.49) 


The new value for in Fig. 3-10) is obtained by calculating the angle a 

from the law of cosines 


a = arc cos\ 


(^t \ t - (Sp sin„)‘ 
2 r, 


and 


and 


(3.50) 


(j)^ = <})^ + a if Q - < z 

<j>^ = <j>^ - a if 9 - > z 


Equations (3.48), (3.49) and (3.50) are used to update the cylindrical 
coordinates (r,h and 0 ) of energy bundles as they undergo scatter or 
absorption/reemission events throughout a single SRB plume. When an 
energy bundle escapes the boundaries of one plume, it is necessary to check 
whether the line of sight of the bundle intersects the neighboring plume. This 
check is provided by representing the neighboring plume as a quadric surface 
in a central coordinate system (Fig. 3- 11). The left plume in Fig . 3 - 1 1 con - 
tains an event site with the line of sight of the energy bundle passing through 
the neighbor plume. The right plume is defined by a quadric expression in 
the central coordinate system. The definition of the line of sight is also 
translated to the central coordinate system and the algebraic expressions 
for the line of sight and the quadric expression for the conical surface of the 
external boundary of the neighbor plume are checked for a possible solution. 

If a real solution exists, there is an intersection of the line of sight with some 
portion of the infinite conical surface. The coordinates of this intersection 
are checked to determine if it lies within the limits of the finite plume. The 
method for representing the plume as a quadric surface is discussed in 
Section 3.5. 
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3.5 SURFACE GEOMETRY PACKAGE 

Heated surface geometries are represented by quadric surfaces. These 
surfaces may be grouped in various combinations to simulate Space Shuttle 
geometries. The external boundaries of the dual plumes are also represented 
as a quadric surface (cone) to facilitate modeling the shading effect of the 
adjacent plumes on the emitted radiation. 

3.5.1 Quadric Surfaces 

The heated surfaces can be described as either quadric surfaces or 
planar surfaces. The surfaces of cylinder, cone, sphere, ellipsoid and 
paraboloid belong to the former type; the surfaces of parallelogram and 
disk belong to the latter. These surfaces can be represented by a general 
quadric equation of the following form 

c, + c, x" + C3 Xj X, + C3 X^ X3 

+ C^X 2 X 3 +C 7 X, tCsXj +C,X 3 +Cjo =0 

The X and X^ are the coordinates of the systenni. The , , , 

are the coefficients that determine the properties of the surface. A detailed 
discussion on these surfaces is found in Ref. 9- 

There are altogether six surfaces, namely: cylinder, frustum, paraboloid, 
ellipsoid, parallelogram and annular disk. On the first three types of surfaces, 
there are optional constraint planes to define the extent of these surfaces. The 
constraint planes are assumed to be perpendicular to the body axis of the sur- 
face. 


The surfaces are defined by eleven parameters in general. These param- 
eters are two radii and the coordinates for three points, as summarized in Table 
3-1, The target surface may be a simple surface as in the cases of ellipsoid, 
parallelogram and disk or a simple surface plus constraint planes perpen- 
dicular to its axis, as in the cases of cylinder, frustum and paraboloid. The 
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Table 3-1 

INPUTS OF TARGET SURFACES 




^2 

^1 

^2 

^3 

C ylinder 

Radius 

N/A 

Center of 
Top 

Center of 
Bottom 

A Point 
not on 
Body 
Axis 

F rustum 

Radius of 
Top 

Radius of 
Bottom 




Paraboloid 

Radius of 
Constraint Plane 
Passing Through 
P2 

1 

N/A 

Vertex 

Center of 
the Base 


Ellipsoid 

N/A 

N/A 

Center 

Pole 

Zero 

Meridian 

Parallelogram 

N/A 

N/A 

Pj, P^, P^ three c 

live corners, counter( 
about the outward nor 

onsecu- 

:lockwise 

mal. 

Disk 

Outer Radius 

Inner 

Radius 

Center 

Arbitrary 
Point on 
Disk 

Normal 

Vector 


N/A; Not applicable. 


3-35 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D496763-I 


basic target surfaces are transformed from the central coordinates to body 
coordinates in which the shape of the quadric surface is readily recognizable. 

3.5.2 Subdivision of Surfaces 

The main surfaces and their constraint planes are further divided into 
sub-areas for more precise heating rate prediction. Four integer numbers 
are used to describe the subdivision. These numbers are NN 1, NN2, NRING 1 
and NRING2 in the program. Their meaning is indicated in Fig. 3- 12. 

The area of each sub area and the position vector of its center are 
computed. Each area is further assigned a ntomber. When a photon tra- 
jectory intercepts a target surface either on the main surface or on the con- 
straint plane, the position vector of the hit point is calculated. Comparisons 
of this hit point position vector to all the position vectors of the sub-areas 
on the given surface are then made in order to identify the proper sub-area 
which was hit. 

3.6 RAY TRACING FROM PLUME TO SURFACE 

The ray tracing scheme can be described as follows. Assume a ray 
Q rig jnating fr om point (N 1,Y1,E1) in the dir ection (U1,U2,U3). Its tra- 
jectory which is a straight line in the three-dimensional space can be de- 
scribed in the parametric form; 

X = XI + U1 . d 

Y = Y1 + U2 . d (3.52) 

Z = Z1 + U3 " d 

where d is a parameter characterizing the distance between an arbitrary 
point (X, Y, Z) on the line to a fixed point (XI, Yl, Zl) on the same line. 

The target surfaces, six in total, can be classified into two categories, 
i.e., the planar surface and the quadric surface. A planar surface can be 
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expressed as 


(X - XO) . nl + (Y - YO) . n2 + (Z - ZO) • n3 = 0 (3.53) 

where (XO, YO, ZO) are given points on the plane; (nl,n2,n3) are the direction 
cosines of its normal and (X, Y, Z) are the general points on the plane. If the 
straight line, as given in Eq. (3.52), intercepts the plane, there will be common 
point (X, Y, Z) such that 

(XI - U1 • d - XO) . nl + (Y1 - U2 * d - YO) . n2 + (Z1 - U3 • d - ZO) . n3 = 0 
therefore 

, ^ (XI - XO) > nl 4 (Y1 - YO) > n2 + (Z 1 - ZO) * n3 
^ “ U1 • nl + U2 . n2 + U3 * n3 ^ 

The position vector of the hit point is obtained by substituting d obtained 
in Eq, (3.54) into Eq, (3.53). A closer look at Eq, (3.54) reveals the following 
obs ervations : 

a. When the denominator vanishes, there is no solution. This is the 
case in which the trajectory is parallel to the plane. 

b. When the numerator vanishes, it means that the given point of the 
trajectory is already a point on the plane. 

c- When the distance, d, is a negative number, it means that the 
direction of the trajectory is in the sense away from the plane 
from the given point. 

Hit points on the quadric surfaces are calculated in a similar way. 

To illustrate the procedure, a cylindrical surface of radius R is considered. 
Let the surface be given in its body coordinates as follows: 

(3.55) 

The common solution of Eq. (3.55) and Eq. (3.52) leads to 
(XI + U1 • d)^ + (Y1 + U2 • d)^ = R^ 


EEPRODTICmiLITY OF THE 
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Thus 

where 


I 




-B+^B^-4*A*C 
^ = 2A 

A = (Ul)^+ (U2)^ = 1 -(U3)^ 

B = 2 (U1 • XI + U2 • Yl) 

C = (Xl)^+ (Yl)^ - 


(3.56) 


If the discriminant, B - 4 AC, is less than zero, there is no solution. This 
means the trajectory does not intercept the target surface. The sign of the 
square root is chosen in such a way that d is the positive and the smaller of 
the two real roots. When d is computed, the hit point on the target surface 
is readily known by substituting d into Eq.(3.52). 

• When there is more than one target surface in the computation, the 
shading effect needs to be considered. The photon ray, being a straight line 
and absorbed by the first intercepting surface, stops at the nearest surface 
it intercepts. The mathematical expressions of all the possible intercepting 
points, however , are calculated in the program. A comparison is made to 
pick the shortest distance and therefore the correct hit point. After the hit 
point is determined, a further check is made to determine the sub-area on 
the surface that is hit. The number of hits of that sub-area is then incre- 
mented. 

The algorithm of the dual plume simulation can be described as follows; 
the two plumes are assumed to be identical in every respect except their loca- 
tion and orientation. Each plume has a local coordinate system which is 
centered at the center of its exit plane with XT axis coinciding with the plume 
axis. In addition there is a central coordinate which is located midway be- 
tween the centers of the exit planes of the two plumes at their zero gimbal 
angles. Each plume can then be located freely with respect to the central 
coordinate system. The coordinate systems and the gimbal angles are shown 
in Fig. 3- 13 and Fig. 3-14, respectively. 
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Let (X 1 , X2, X 3) and (U 1 , U2, U3) be the position vector and the trajectory- 
direction cosines in the central coordinate system. The primed quantities 
(X 1 X2', X3') and (U1 U2', U3') are the corresponding quantities in the plume 
local coordinates. (PP(1), PP(2), PP(3)), and a are the location and gimbal 
angles of either plume. Transformation from the plume coordinates to the 
central coordinate can be done by following matrix multiplications. Note that 
ip is either or ip^ plus k/Z and a is either or as the case may be. 



/cosCT 0 sina 

sinCT sin^i cosip -cosdr simp] 
-sina cosvp sina cosa cosipy 


/Xl'\ /PP(1)\ 

f X2' | + (PP(2) I 
\X3'/ \PP(3)/ 


(3.57) 


U1 

U2 

U3 


( cosa 0 

sina simp cosip 
-sina costp sinip 


sina 
-cosa simp 
cosa cosip 



(3.58) 


Transformations from the central coordinate system to either plume 
coordinate system are the inverse of the above and are given as follows; 




sina simp 
cos 4^ 

-cosa sin>p 


-sina cosipX ^ " PP(1)\ 
simp I ( ' PP(2) j 

cosa cosip / \x3 - PP(3)/ 


(3.59) 




sina sinip 
cosip 

-cosa simp 




(3.60) 


When an energy bundle escapes the boundary of a plume, the location 
of the point on the plume boundary and the direction of the bundle can be re- 
corded on a data tape. When a large number of sample energy bundles are 
collected, typically 100,000 samples, the data will give a good definition of 
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the plume's radiation field. This process is the IPLUME = 3 option in the com- 
puter code. The energy bundle can be traced directly as it leaves a plume. A 
random number is first called to determine in which plume this energy bundle 
is originated. The location and direction are then transformed to the central 
coordinate system. A check must be made at this point to determine if this 
energy bundle will intercept the other plume. If it does, the interception on the 
other plume is treated as a reemission point in that plume, and the tracing process 
described in preceding sections is continued in that plume until the energy bundle 
again leaves the plume. If the energy bundle leaves a plume and does not inter- 
cept the other plume, then the geometry package is called to determine if it hits 
any part of its surfaces. This procedure is the IPLUME = 1 option in the code. 

If the data tape is used as in the case of IPLUME = 4 option, instead of 
generating the energy bundles within the plumes, the presence of the other 
plume is treated as a shading surface. When the energy bundle hits the bound- 
ary of the other plume, the energy bundle stops there. Otherwise a check is 
nnade to determine the hit point, kloreover, due to the assumption that the dual 
plumes are identical when the location and direction of an energy bundle are 
read from the data tape, the same set of data applies to both plumes. There- 
fore when a set of data is read from the tape, two rays which are originated 
at each of the two plumes in a mirror symmetrical position are calculated. 

In the case when the plumes are gimbaled, this symmetry still applies, be- 
cause the location and direction are given in terms of the plume local coordi- 
nates which are not affected by the gimbal angles. 

3.7 HEAT FLUX EXPRESSIONS 

When an energy bundle escapes the boundaries of the dual SRB plumes, 
it may or may not strike a surface of the geometric model of Space Shuttle 
surfaces. The methodology for discriminating whether the bundle intersects 
a spacecraft surface has been discussed in Section 3.5. If the bundle hits a 
surface, the energy contained in the bundle is absorbed by the surface. The 
radiant heat flux to a given surface is ultimately determined by tabulating the 
number of energy bundles that strike the given surface relative to the total 
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number of bundles used in calculating the distribution of thermal radiation 
leaving the plumes. 

The amount of energy contained in an energy bundle for a specific SRB 
plume is dependent on the number of samples, N^, and the number of absorption/ 
reemission events that occur during the simulation. The total number of samples 
plus reemissions multiplied by the energy per bundle must equal the total energy 
emitted by all regions of the plume as given by Eq. (3.11), as 

(N + N ) e 
' e re 

Equation (3.61) defines the energy per bundle for a specified sample size, N^. 

The heat flux to a given surface is given by 

qn = e/Ag (3.62) 

where is the area of the sub-surface and N, is the number of hits in that 


(N CT eJ 

+ K E^ 

1 a bl 

a b 

1 


V. 


(3.61) 
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4. VIEW FACTOR MODEL 


4.1 HEAT TRANSFER CALCULATIONS USING VIEW FACTORS 


The radiative heat exchange between surfaces can be calculated with 
the help of the view factor concepts. The view factor, hF^ which is also 
known as configuration factor, shape factor and other names, is defined as 
the fraction of energy leaving black surface dA^ that arrives at black surface 
dA^ (Fig. 4-1), 

cosp. cosp^ 

dFi_2 = dA^ (4.1) 

rr S 


Equation (4.1) shows that the view’ factor depends only on the size of 
dA^ and its orientation with respect to dA^. The differential notation dF indi- 
cates that the view factor is for energy transfer to a differential element, dA^, 
from a differential element, dA^, The view factor for two finite areas, A^ 
and A^, is calculated by integrating over both A^ and A^ 


F 


1-2 



If 

AiA^ 


cos(3^ cosp^ 


7T S" 


dA^ dA^ 


(4.2) 


The view factor F^ j which is defined as the fraction of energy leaving 
black surface that arrives at black surface A^, can be given similar to 
Eq. (4.2) 


2-1 


kff 

A2A1 


cosp^ cosp J 

5 T s 


dAj dA^ 


(4.3) 
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Therefore, the reciprocity 

relation 

follows 


AfF 

1-2 "" ^ 

F 

2 2-1 

(4.4) 

The radiation exchange from A^ to 

A^, ^2’ from A^ to 

Af, 

q , can be given immediately. 




^2 

= a A^ 

4 

F T 

1-2 M 

(4.5) 

*11 

= crA2 

F t"^ 

^2-1 2 

(4.6) 


The view factors, and depend on not only the area and orien- 

tation of the surface but also on the shape and local curvature of the surface. 

The view factors of some standard configurations are found in the literature. 

For those configurations for which the view factors are not available, the 
process to compute the view factor by evaluating the integration given in 
Eq. (4.2) is difficult in general. 

4.2 MONTE CARLO METHOD 

Since the view factor is defined as the fraction of total energy leaving 
one surface that arrives at the other, it is possible to divide this total energy 
into many bundles which are emitted diffusely from the emitting surface. The 
fraction of the total bundles that is intercepted by the other surface would be the 
view factor F^ 2 ' this manner the view factor is obtained without carrying 
out the lengthy integration as given in Eq. (4.2). This is the underlying idea 
of the Monte Carlo method. 

The Monte Carlo method is especially powerful when the geometry of 
the exchange surfaces is complex. There are, however, assumptions in- 
volved. The assumptions required when using the view factor to compute 
radiant heat exchange are that the surfaces involved are diffuse-gray emitters 
and reflectors, that each surface is isothermal and that the total flux arriving 
at and leaving each surface is evenly distributed across the surface. The 
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assumption required with the Monte Carlo method is that the reflector surface 
is of some reasonable finite size so that the number of hits recorded on the 
surface would be statistically meaningful. This assumption is not as restric- 
tive as it appears. In the case where the reflector is small the view factor 
can be calculated with reciprocity relation by first computing the reverse 
view factor from the reflector to the emitter. However, there is no apparent 
criterion to distinguish a small target surface from a non -small target sur- 
face. It is a matter of how large a solid angle that the target surface subtends 
and how many emitting bundles are used. In the present code, the view factor 
and its corresponding reciprocal are computed at the same time. These are 
given in the output side by side for convenient reference. 

4.3 DISTRIBUTION FUNCTION 

In the Monte Carlo method the energy bundles are emitted from a sur- 
face uniformly in a cosine distribution. The process involves picking an 
emission site on the surface and choosing a direction of the trajectory. When 
a large number of energy bundles are dispensed, the emission sites will be 
distributed uniformly over the emitter surface and the trajectory directions 
will be distributed according to the cosine law. The spatial distribution will 
be discussed in the next section. The directional distribution is discussed 
here. 


The direction of the trajectory can be determined by the angles r) and 9 
(Fig. 4-2) when the location and the local normal are given on the emitting 
surface. When a large number of emissions are chosen, the trajectories will 
be distributed according to the cosine law about the local normal. 

Let N be the total number emitted from the source point and dn the 
e 

number emitting through dA. Since the number density is proportional to 
cosr], the number of trajectories lying within du is 

du = sin r; d T) d0 
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d N = ^ N COST] 

^ e 4jr m ' 


N 


= sinr) cos r) dr) d9 


where N is the maximum number density. The total number 
m 


N 


t/ 2 Ztt 

- I I sinr) cosr) d9 dr) 

0 0 


N 


m 


The joint probability density function is therefore. 


^ = fd9 dr) = — sinr) cosr) d0 dp 
N ' r 


The marginal distribution functions B^nd H ® 


respectively are 


7t/2 


( 0 ) 


= y* ^ sinp cosp dp = -^ 


0 

Ztt 


(n) 


/ i =‘"'’ 


cosp d0 = 2 sinp cosp 


0 


Therefore the corresponding random numbers are 

^9 =/ ^ 


u 


r] ^ 

= y* 2 sinp cosp dp = sin p 
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which moans 

9 = 2jr • Uq 


(4.7) 


r] 



(4.8) 


Equations (4.7) and (4.8) are the relations between the random numbers 
Uq and drawn from the random number generator and the values of 9 and rj 
used to determine the trajectory. 

4.4 EMITTING SURFACE GEOMETRY PACKAGE 


The geometry package in the view factor calculation involves two parts. 
The first part is the geometry of the emitting surface which will be discussed 
here. The second part is the geometry of the target surfaces which is the 
same package used in heating rate computation and was discussed in the pre- 
ceding section. 


The emitting geometry package includes nine surfaces, i.e., half cylinder, 
half frustum, hemisphere, parallelogram, annular disk, cylinder, frustum and 
sphere. The input and output of these surfaces are all done in a subroutine 
VFEMIT, When the subroutine is called to output a trajectory, a point on the 
emitting surface is first chosen as follows: 


1 . Cylinder 

HELV = H • RNl 
TH = (0.5 - RN2) . t/2 


where 

HELV is the height above the base 
H is the height of the cylinder 
TH is the azimuthal angle, and 
RNl and RN2 are two random numbers. 
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2. F rustum 


HELV = H • (R2 - RL)/(R2 - Rl) 
TH = (0.5 - RNl) . jr/2 


where 

Rl and R2 are the- radii of the base and the top of the frustum, 
respectively, and 

RL = R2 * ^RN2 . 

If RL < Rl, another RN2 is drawn. 


3. Sphere 

A random direction is chosen. The point on the spherical surface is 
determined by a trajectory starting from the center of the sphere in the 
chosen direction to the point on the surface. A final check is made to ensure 
the point falls within the range of the hemisphere. If the point on the spherical 
surface is not within the given hemisphere, another random direction is gen- 
erated and the process is repeated. 

4. Parallelogram 

D1 = RNl • H12 
D2 = RN2 • H23 

X = PI + C 12 • D1 + C23 . D2 

where 

X is the position vector of the random point on the plane 
PI is the position vector of corner PI 

H12, H23 are the distance PI to P2 and P2 to P3, respectively, and 

C12,C23 are the direction vectors of the lines PI to P2 and P2 
to P3, respectively. 


4-8 


LOCKHEED • HUNTSVILLE RESEARCH A ENGINEERING CENTER 



LMSC-HREC TR D4 96763-1 


5 . Annular Disk 

TH = RN2 • Ztt 
RR = R2 .^RNl 

where 

TH is the azimuthal angle subtended by the random point, and 

RR is the random radius. If RR is less than the inner 
radius Rl, then another RNl is generated 

Now that the random point on the emitting surface is chosen, the direction 
of the emission is the next to be determined which is described in the following 
paragraphs. 

4.5 RAY TRACING 

When a random point is chosen on the emitting surface, the direction 
cosines of the local normal at that point are known. The direction of the 
emitting trajectory is then generated by random numbers in such a way that 
the direction would fall into a cosine distribution about the local normal. A 
set of direction cosines (UT, U2^ U3') with respect to the local coordinates 
are generated as follows. 

9' = • RNl 

4>' = sin"^ ^RN2 

Ul' = cos^' 

U2' = sin^' COS0' 

U3' = sin^' sin9' 

The local coordinates which coincide with, the local normal at the chosen 
surface point can be related back to the central coordinates. The direction 
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cosines (U1 U2’, U3’) are transformed into (Ul, U2, U3) which are based on 
the central coordinates. Let the local normal be (N1,N2,N3), The trans- 
formation is: 


U1 


cos^ o sin^ 


ur 

U2 


sin^ sin0 cos9 -cos^ sin9 


U2' 

U3 


-sin^ COS0 sin0 cos^ cos9 


U3’ 


where 

^ = ACOS(Nl) 

9 = 7t/2 + ATAN (N3,N2) 

Now that the location of the starting point and the direction of emission 
are known, this information can be fed into SORTNG subroutine to determine 
if and where the target surface is hit. The number of hits are recorded 
according to the sub-areas on each target surface. 

When the specified number of trajectories are dispensed, the view factor 
and its reciprocal view factor of each sub-area are computed b^ dividing the 
number of hits by the total emissions. The results are grouped together by 
each target surface in the output. 
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5. EXAMPLES AND RESULTS 

5.1 RESULTS CONCERNING PLUME CHARACTERISTICS 

The two plumes from the two SRBs are assumed to be axisymmetric 
and identical to each other. The radiative power of the plumes varies with 
many parameters, e.g., trajectory point, particle size and density, after- 
burning, gas constituents, etc. The influence of these plume flowfield prop- 
erties is reflected in the Monte Carlo solution to the radiation transport 
problem. 

The energy bundles are traced within the plume until the energy bundles 
leave the plume boundary. Data tapes are generated recording the direction 
and the location of the line of sight of each energy bundle leaving the plume. 
These data tapes are used to define the line of sight of energy bundles leaving 
the plume in the radiation heating rate calculation. This circumvents the 
necessity of performing the Monte Carlo simulation for every computer run 
and provides a substantial saving in computer execution time with no sacrifice 
of sample size or accuracy. 

Distributions of the energy bundles as they leave the boundary of a single 
plume are plotted on the following pages. Figure 5-1 shows the distribution in 
the Xj direction which is along the axis of a single plume. Figure 5-2 shows 
the distribution of the bundles leaving a single plume relative to the angle r) • 
The rj angle is the angle between the trajectory of the energy bundle and the 
Xj axis. Distribution of energy bundles in the (j> direction is always uniform. 

5.2 EXAMPLES OF HEATING RATE CALCULATIONS 

The radiant heat flux to a disk at a distance -0.1 over a single plume 

for sea level and 72,000 ft trajectory point s is plotted in Fig. 5-3. Figure 5-4 
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shows similar data for a single plume except the disk is -4.94 R above the 
exit plane, which is about the same elevation of the bottom of the ET. 

The effects of gimbal angles and different altitudes on the radiant heat 
fluxes are shown in Figs. 5-5 through 5-10 for the dual SRB plumes. The 
heating rates are shown for a rectangular plane extending from = -6 to 6 
and = -5 to 5 over the dual plume. Figure 5-5 shows the heating rates 
along the strips parallel to X^ axis. The two obvious dips along the center 
strip are due to the shading of the nozzles. The distance between the rec- 
tangular plane and the nozzle exit plane is denoted by X^^. Figure 5-6 shows 
the similar case, except the plane is located at X^ = -4.94 where the plane is 
tangent to the bottom of the external tank. The heat flux is more uniform at 
X = 4.94 and the intensity is reduced by a factor of 2 to 3 as compared to the 
X = 0 location . 

Figures 5-7 and 5-8 show the effect of gimbal angle for the sea level 
plumes. Figure 5-7 shows the heat flux on strips of the rectangular plane 
when the plumes are tilting 30 deg toward the X^ axis, i.e., a = -30 deg, 

= 0 deg. The heat flux on the side where X^ > 0 is a marked increase over 
the side where X^ < 0. The heat flux on the strip toward which the plumes 
are tilting is 3 to 10 times higher than that on the strip away from the 
plumes. Figure 5-8 shows the heat flux on strips of the rectangular plane 
when the plumes are tilting 25 deg toward the X^ axis, i.e., c = 25 deg, 4* = 90 
deg. The heat fluxes are lopsided to the same direction the plumes are tilting. 

The heat fluxes along one strip near the edge of the rectangular plane 
are plotted in Fig. 5-9 for the 9430 ft plumes. The four curves represent the 
following cases; (1) a = 0° 4 ^ = of X^ = 0; (2) cr = 0? i]; = 0? X^ = -4.94; (3) a = 25? 
4 = 0? X j =0; and (4) a = -30? 4 = 0? = 0- It is seen that the gimbaling of 

the <j angle would increase or decrease the heat flux at this particular location 
by a factor of 2. 

The heat fluxes for no gimbal for the 72,000 ft plume are plotted in 
Fig. 5-10. The bottom curve is the heat flux at X^ = -4.94. The heat flux 
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Fig. 5-7 - Heating Rate of a Rectangular Plane 
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at Xj = -4,94 is seen to be only one third of that at X ^ =0 for the same row 
of subtargets. 

A simple geometry package (see Fig. 5-11) was put together to compute 
the heat flux at the end dome of the ET. The shading effects of the nozzle are 
simulated with a disk placed at the nozzle exit plane. The end of the ET is 
divided into three rings which are further divided into four parts. The heat 
flux on each division is summarized in Table 5-1. Table 5-2 shows the heat 
flux when the sea level plumes are gimbaled in pitch plane and in yaw plane. 
The effects of gimbaling are especially pronounced when the plumes have large 
pitch angles (a large, ijj = 0). Table 5-3 shows the heat flux of gimbaled dual 
plumes at 72,000 ft altitude. It is observed that the gimbal angles in both the 
pitch plane or yaw plane increase the heat flux at the ET dome substantially. 

Two models were made to facilitate the study of heat flux within the 
nozzle shroud enclosure. One model is shown in Fig. 5-12. Four target sur- 
faces constitute the model. Surface 2 is the skirt of the nozzle plus the end 
constraint plane and a stiffener. Surface 4 is another stiffener at the mid- 
height of the skirt. The heat fluxes on these two surfaces are given in Table 
5-4 for a single plume at three altitudes. 

Another nozzle enclosure model is shown in Fig. 5-13. This model is 
in closer agreement in its dimensions to the actual shroud design on the SRB. 
The target surfaces are designated by numbers. The heat fluxes are given 
in Table 5-5 for a single plume at three altitudes. 

The geometry model of the ET and SRB structure is shown in Fig. 5- 14. 
The target surfaces are designated by numbers. The end dome of the ET is 
divided into four rings on its surface. The plumes at three altitudes are not 
gimbaled. The heat fluxes are summarized in Table 5-6. 
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Table 5-2 

EFFECTS OF GIMBAL ANGLES ON HEAT FLUXES 
AT ET DOME FROM DUAL SEA LEVEL PLUMES 


(deg) 

0 

90 

0 

90 

0 

o (deg) 

0 

8 

8 

35 

35 

1 

1.14* 

1.13 

0.32 

1.78 

0.81 

2 

1.14 

1.62 

0.32 

1.62 

0.81 

3 

1.46 

0.97 

0.97 

1.46 

2.27 

4 

1.46 

1.14 

0.97 

2.27 

2.27 

Avg. 

1.30 

1.22 

0.65 

1.78 

1.54 

5 

0.93 

1.12 

1.31 

1.68 

1.31 

6 

0.93 

1.87 

1.31 

2.24 

1.31 

7 

2.43 

1.12 

1.31 

2,06 

2.62 

8 

2.43 

0.75 

1.31 

2.62 

2.62 

Avg. 

1.68 

1.22 

1.31 

2.15 

1.96 

9 

0.91 

1.83 

2.05 

4.34 

2.05 

10 

0.91 

1.83 

2.05 

2.51 

2.05 

11 

1.83 

3.65 

2.97 

4.11 

4.79 

12 

1.83 

1.83 

2.97 

4.11 

4.79 

Avg. 

1.37 

2.28 

2.51 

3.77 

3.42 


* 2 
Heat flvix in Btu/ft -sec 
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Table 5-3 

HEAT FLUX ON THE ET DOME WITH GIMBALED PLUMES AT 72,000 FT 


Targ et 

o (deg) 

0 

35 

35 

(deg) 

0 

vU __ 

90 

0 


1 

0.67’ 

0.88 

0.60 


2 

0.67 

1.00 

0.60 


3 

0.70 

0.88 

1.33 


4 

0.70 

1.05 

1.33 

Avg. 

0.68 

0.96 

0,97 


5 

0.91 

1.57 

1.01 


6 

0.91 

1.54 

1.01 


7 

1.03 

1.42 

1.80 


8 

1.03 

1.36 

1.80 

Avg. 

0.97 

1.47 

1.45 


9 

1.43 

2.03 

1.78 


10 

1.43 

2.37 

1.78 


1 1 

1.28 

2.32 

2.69 


12 

1.28 

2,37 

2.69 

Avg. 

1.36 

2.27 

2.24 


* 


Heat flux in Btu 
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Table 5-5 

HEAT FLUX ON A NOZZLE SHROUD MODEL 
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Fig. 5-14 - Model of the External Tank and Solid Rocket Boosters 
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5.3 EXAMPLES OF VIEW FACTOR CALCULATIONS 

Several examples of view factor calculations are presented in this 
section. The first example involves the geometry of three disks, as shown 
in Fig. 5-15. The surface 1 is the emitting surface. 

To compare the results with the exact solution and the results obtained 
with different sample sizes, the view factor of the entire surfaces is used 
although the view factors from the emitting surface to the sub-areas on the 
target surface are also available. The exact solution of view factor ^-^_2 
given as follows (Ref. 6): 


R^ = r^/h, R^ = r^/h 

X = 1 + (1 + 

Fi_2 = ^ ^ ^ {R^/R^)^)/Z 

The results are summarized in Table 5-7. 


A second example is given in Fig. 5-16, where a hemispherical emitting 
surface of unit radius is facing a dish-shaped target surface. The exact solu- 
tion of the view factor calculation is available in the literature (Ref. 6) and is 
given as follows: 



Comparisons of results are given in Table 5-8. 

The third example demonstrates the heat flux calculation to a point in 
the base plane and a point on the base of the ET using the view factor method. 
A small disk with a radius of 0.001 is placed at either position 1 or position 2 
at the end point of the ET (Fig. 5-17). The two SRB plumes for the 72,000 ft 
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Fig.5-15 - First Example of View 
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Table 5-7 


COMPARISON OF RESULTS FOR EXAMPLE 1 


View 

^\Factor 





Sample 

Method 

F 

1-2 

^2-1 

^1-3 

^3-1 

1000 

0.2840 

0.0710 

0.0370 

0.00925 

2000 

0.3020 

0.0755 

0.0345 

0.00863 

5000 ' 

0.2950 

0.0738 

0.0322 ■ 

0.00805 

10,000 

0.2871 

0.0718 

0.0336 

0.0084 

20,000 

0.291 1 

0.0728 

0.0358 

0.00894 

Analytical 

0.2918 

0.07295 

— 
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Table 5-8 

COMPARISON OF RESULTS FOR EXAMPLE 2 


View 

^'v^actor 

Sample 

Method 

r 

h 

^1-2 

^2-1 

Analytical 

1 

5 

0.01942 

0.03884 

7500 

1 

5 

0,01960 

0.03920 

Analytical 

5 

5 

0.29289 

0.02343 

7500 

5 

5 

0.2944 

0.02355 

5000* 

5 

5 

0.2896 

0.09267 


5 ^ 

Radius of emitting hemisphere = 2. 
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trajectoTy point are represented by two frustum bodies, Tbe average emis- 
sive power of the plume is assumed to be either 40 Btu/ft -sec or 25 Btu/ft 
sec. The results are summarized in Table 5-9. The schematic of the SRB 
plumes and target points is shown in Fig. 5-17. 

Asa demonstration case, view factors from a sphere and a cylindrical 
rod (simulating a strut) to the surfaces within the nozzle shroud enclosure 
shown in Fig. 5-12 were computed. The orientations of the sphere and rod 
relative to the shroud are shown in Fig. 5-18. The surfaces are designated 
as target 1, target 2 and so on and their orientation is seen in Fig. 5-12. The 
ej^itting surfaces are the sphere and rod and are designated as number 1 and 
2, respectively. The computer output is reproduced in Table 5-10. To limit 
the data printout to a reasonable number of pages, only data for target 1 (the 
nozzle) and target 2 (the conical shroud) are shown. 
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Table 5-9 

HEATING RATES BY VIEW FACTORS FOR 72,000 FT PLUME 


Location 

1 

. . . . 

2 

View Factor ^ 

View Factor F^ ^ 
Sum F - 1 

1.9329-09 

1.8959-09 

3.8288-09 

6.0481-10 

6.0736-10 

1.2122-09 



2 

Heat Flux, Btu/ft -sec 

at 40 Btu/ft^-sec 
2 

at 25 Btu/ ft -sec 

28.28 

17.66 

8.95 

5.59 

Emitting Surface: r = 0.001 ; A ~ 3.14159 x 10 ^ = 1.1575 x 10 ft 

* ^2 

Surface Area 2 or 3 = 580 = 21370 ft 

Total Emission from 2 or 3 = = 21370 x 40 = 8.55 x 10^ Btu/sec 

E^ = 21370 X 25 = 5.34 x 10^ Btu/sec 

Heat Flux at 40 Btu/ft^-sec = 8.55 x 10^ x F/A 

Heat Flux at 25 Btu/ft -sec = 5.34 x 10 x F/A^ 


Non - dim en s ional . 
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Table 5-10 

VIEW FACTORS FOR SHROUD GEOMETRY 


input data roH target nu. I 


7 I ? 

-1 


J I I 

- f r <4 * ' p 


.ooro 

0 c! r n J 0 • 0’6 do 


*. H soar 


^ ?: 'i r 0 p h A T rGTT-M JrTTfTT 

. I *1 1 • " - r : 3 

;• ''3 • \z^>^TZ*r\ 

.p*;: ' -<c 


• r : 

'•re ^i^nc 
« KCO-'O^Gt 


f ^ M r f M. r K A c OUE 
1 'y - . i H ♦ '"T 
- i ■ ’ ’ - ♦ C I 

. I c ^ ‘ ‘ ♦ C I 
• C ^ ^ 


c. .• ^ 

C a 
c H U 

C*. :« 
C 6 ■ 

C 7 - 

C « ' -« 

c - - 


/ ^ I X J w j 


c j . . t . I : " ^ I 

i< f‘ i ‘u» ui ^ A I A I . I-! n ouD u <‘Nf) A r'A y . i c ?Y+ oT T^fScUTTc ♦d(T 


ui 

I 

OJ 


CU} I 

c I I 
Cl 7 
C I T 
cf*i 
C A 


> IC It.-.r 
’ • I '• 


r':4'“ 


E. f OH cnrjSTRAiNT planes 

r - ^ I - -- *- 




y 1 '^cuc t 

M \ K M T S 

1 • 

*■■’0 1 NTS 

FOP' VT 
VP 1 » 

ON constraint 

DISK T( ,'IOSOC 

•OOCOO I.OOOOC 

> S( 

•00000 1,00000 ,00000 ) 

'Z-1^ 'POlfTT' V'FCT^ 
-.7S-S^^4CC 





tRS 

VP7- .23Sl VV4.0? 

VPJ* •a773834.Q0 

AREA- 

.2532B3^0o 

— 77 “ 

“v PT»“ 

- . ? s : : ^ ^ : j 

7p77 • 6 4229'^ ♦0? 

.6H22894-n3 

AREA- 

. 2T32"53*oo 

3 • 

VIM • 


VPZ* •B773B3+07 

VPJ« •23Sr94>3D 

area- 

. 2S32834-00 

M . 

V I * 1 * 

- • 7 b ^ ♦ ,M ’ 

VP2» •B77JB3<»*aO 

'vP3» -.23S09R4*C10 

AREA- 

.253283*00 

S . 

V TM • 

7S - r ' 

VP2a •6V2?y 

VP3* ‘••692289^T3 

AREA- 

.253283*00 

6 • 

V IM =* 

-.?s 

‘ V p 2« “ 'V2 3s : V V ♦ Or 

'‘VP3* -•8773834nO 

AREA- 

• 253283-»-0D 

7 . 

vfM ■ 

• . ? s : ? r j ♦ } r: 

VP2a -.23S^9q4.or; 

VP3« ••977383>3T 

AREA* 

.253283*00 

' ■ n V 

V p 

5“ 

VP2* -•6M7?89>00 

~vT>3i '-.6«2?09*OO “ 

ARE A - 

•iS3 283+00 — 

9 • 

VP 1 « 


VP?« -.B773B3+00 

VP3* -•23S:994'T3 

area- 

• 2^32«34'00 

10. 

VI' 1 * 

• /s"7- iVno ■“ 

“^vT> 2 * “-T;T7 7 3 B'3 V 

”VP3» •2 3So9 4-*-'>0 " 

A R E A - 

• 2^32d3-»-i3D 

1 1 . 

V fM * 


VP2« -•6922694.00 

VPJ» •69228V4-n0 

ARE A- 

.253283>0c 

17* 

V FM « 


" vp?» ■-.?3s:99*ao ■ ^ 

VP3* •8773e34'lt) 

area-" 

. 253283^00 

1 3 • 

VP 1 * 


VP2« • 1 B76HM4Q0 

VPJ- •7?0296^:»0 

area- 

.202161 >00 

” ‘ 1 ti . 

v ’P T ^ 


’VP?- 

"^P3-‘*''.5»26S2>3^ “ 

A FTE A o 

. 2'’2 1 S T*0{5 

1 * 

V P 1 » 

- . 7 S " 7 ♦ .1 T 

VP2» • ;r.:296*00 

VP3* *|87694^3'? 

AREA■ 

.202161*30 

« '> • 

vn* 

- . 7 s ^ ^ ♦ ; - ■ ' 

vP 2= ~»7‘rT2 9A + dd 

“vP3» '-• 1 B7699 + *n 

" ‘ A«t A- ' 

.202 1 6 1 *Oo 

i 7 . 

V P 1 = 

-- 7^- “ : /. 

V P 2 « • S 1 7 6 S ? ♦ 0 J 

VPi* -tSl 26S2^ TO 

A«EA- 

.202161*00 

1 B« 

\l p 1 » 

;s- 

V P 2 • 1 B 7 6 9 M ♦ J nf ' 

vP -• 7''0296^-1? 

AREA- 

. 2 C 2 1 6 1 *Oo 

IV* 

^ P 1 a 

- . / s : 1 ♦ 

VP?» -•IR7694*Q0 

VP3« -•7^r?296^'O0 

aRf.a- 

.202161*00 

' ■ 7u.' ’ 

-vP 1 » 

-*7s - j: . *.: : ■■■' 

V P 2 S 1 2 ^ ^ 7 ♦ 3 0 

■'VP3« -•Sl2652>J*' 

aRea- 

.2''2l6l*03 

> 1 . 

. P 1 = 

- . 7 "1 1 * '* 

V p 7 - • 7 '' 2 V (1 ♦ C 

VP3» 1 8769H4-0? 

ARtA- 

.202161*00 

7 . 

- * t = 


- • 7^' ?'M, *0S 

’vPj« 

A H L A * 

. 2?2l6l*3D * 
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23 . 

VP 1 • 


VP2" 

-»5| ?65?*0e 

VP3« 

•Si 26S2^00 

AREA- 

.2021.1*00 

24 . 

VP 1 « 


Vp2« 

-. 1 07644-1.00 

VP3« 

•700296*00 

AftCAa 

,2FFiTToc 

2S . 

VP 1 ■ 


VP 2« 

, I 4 1 94 + 00 

VPJ- 

•5?32l n+TT 

AREA- 

• |Sl0H0*0Q 

2a. 

VH 1 « 

-• t2sr^'::*r;i 

VP 2 « 

. 303D 1 6^03 

VP3* 

•30301 6*00 

AREA- 

* ISIOHO^OO 

27 . 

VP J ■ 

1 

Vp2* 

•S2321:^0n 

VP3* 

• 1 4CI 944..30 

AREA- 

. 1 0‘<0*00 

2 H • 

V P 1 * 


VP 2» 

.b2l2l:*00 

VP3s 

-• 1 40194*00 

area- 

. 1 s 1 o.*r*?G 

79. 

VP 1 * 

1 2^1'^^*'' 1 

VP2“ 

.393:J64-00 

VP3s 

-•3ft3ci6*00 

AREA- 

* |S|Q4C^00 

30 . 

VP 1 ■ 

- . 1 2 ^ I 

VP2?* 

• 140 1 94^00 

VP Js 

-.S732 1 r*00 

area- 

• iSiOHO^'io 

3 I • 

;p 1 » 

- ♦ I 7 ^ ^ 1 

VP 2a 

-•14TIV44.00 

VP3« 

-•S232I 0*00 

AREA- 

• 1 s 104C^O;^ 

)2. 

V P 1 • 

- . 1 2 ‘.1 0 - J ^ 0 I 

VP 2 a 

^ 3a 30T6^CO 

0P3« 

-•303016*00 

area- 

,|SlO40*0c 

33 . 

j/ p 1 «= 


VP 2s 

-•92321? >00 

VPja 

SlZ > ^3 1 9 4 •*• 0 0 

area- 

• IS10*0*2C 

5 4 

V P 1 * 

r V ♦ ^ 1 

V P 2 a 

-.b?32i:^00 

V P 3 * 

• 1 4:? i 9 4*00 

AREA- 


5f>. 

V J’ 1 « 

1 2'^^^ j^:ii _ 

yr.2a 

3430 1 6>0*? 

yP3" 

•303016*00 

AREA- 

_ * »040^2C 

3a. 

V p 1 ■ 

I 2bo:*^^ * 

VP2s 

- • 5 4P I 94*00 

VP js 

•S232IC*30 

AREA- 

.ISi040*Ol” 

IINS j’h A r 

I'r Oi l 

k" no .1 








K M f : A ' 1 

A PE A 







i 3 ’ * S 3 1 ! H 4 - - I 


U rUNSthA!’*T DISK 'no*? “ 

. •< Hr A N ,A A 

I , .V6i7P9 + or 
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Input data for target nu, 2 

2 12 8 9 7 3.C'0r?^ .croo 2 .S 0 CH 

-7,SD0P ■ rocrr^ 

T r: JV a TT'U’J 'fTA'TTrx ~ 

• 1 I • D t* J " T • ;? C 3 0 n 

. 1 -:* '' ' .' 1 : :: : 'j^n't rrco“n >'0 

.icodgo^oi 


2 . 2 St ?3 

■ZV^CTTJ “ ZT'ZSOa 


V. Ok'f- M C 1 l r S FOR A CONE 

c r*= ~ -vi 4 

C? ■ .1 • ^ _ 

c 3 » • i : ^ 

Cl 3 

c <s = 

C 7 i 

c H * 3 

c V 3 . c T :; c» 

■ •' ;G r * - • I G I ^02^r2 



it ( Mi A 1 

A > • RDu.JtJU INO A T A f 


♦Cl .22SOCC>Ol 
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C’>r.^F IC It!iTS FOR COUSTRAiNT PLANES 







CM ^ • 1 

C 1 ’ • 



» -1 - -V ■ 








C 1 » * 

c t. ’ i I • ■ 










c '1 i ■ 1 r ■ * - p 
itf L f i R t . N C L 

(- ^ -j 

k)INTS 0'^ constraint 

DISK 

Ti -.5030C 

•00300 3.0O00D I 

S( 

•SOOOO 3.00000 

• 00000 1 












CO'ioonr^jTs 

1 • 

r 0 R Ti r’Y p 0 T NT' vcrro rs 
vPi* -.o2g:’‘^ + :3 vp2« 

•76M32S^0D 

VP3- 

• 2s&zsa'»oi 

AREA- 

.4I2Q&0+00 


O ^ 

3 » 

VP ! « - p/>2SG ".>00 

7p2^a“ 

VP2« 

‘ .znayrrvri 

. 2RS750 ♦O i 

"vn» 

VPJ* 

“TjAsa'i 7*3 r 
.7AH3Z5^m 

AREA* 
ARE A ■ 

• 4 1 28^0*00 
,4l28S0^00 


go 

'pa 

H • 

s# 

V P 1 ■ ‘ - >1 2 G 0 " w ♦ G 3 
VPI* -.*2S.?'^J^03 

p 7 ■ 

-"■.‘2nV2^>0T 

.208817+01 

"VT» j a “ 

VP3« 

■-.7iH32S*nO 
-•20081 7+01 

aRl A« 
AREA- 

• 4 1 afisc-^oc 

.41 28S0^00 


6 • 
? . 

V p 1 = '- .62G::G''*Ga“ “ 
V P 1 ■ - . 2 G 0 > 3 ? 

VP2* 

VP?a 

'■ V7"5m 32^ + 03 
-.7Aq32S+00 

'VP3»' 

VP3« 

-•28525.C*01 ■ 
-.28S2S3*01 

area-"* 

AREA- 

.4l28bC>30 
,41 28&r» + 00 


O 

3 

rt * 

9 « 

— V 1 •■”".■«■? 5 ; ? : ♦ c 6 
vt‘i» 

VP?» 
VP 2* 

• . 2O8TIT+0 1 
-.2nS2S'" + Oi 

VPJ." 
vP J« 

-•2^flQ I /♦Ot 
-•7AR32S+30 

area- 

area- 

r285G>Oo 

, 4 1 28&0+00 


1C. “ 
t ) • 

V n 1 a - i ii I* S 0 + G ^ 

VPl* -.A2ST'^^v''3 

‘”~”vp 2*^ 
vp2 ■ 

-. 2''rt81 7 + 01 

“vTTJi" 

VPO- 

“ *7iR32S^00 
.70881 7>3l 

ARE A- 
ARE A- 

im 785t?*P»o 
. 4 1 2B50*0C 



12. 
1 3. 

-"’VP I . i2G0:j*GO 

V p 1 a - . H 7 s 0 ci 

V P 2 ■ ’ 

VP 7 a 

“*7^m32S+6o 
• 7 M r : 6 ! ♦ C 3 

VP3«' 

VP3» 

■ •28S25>c*D 1 
•?74(9H+nl 

area-” 

area- 

.4 1 285C^0o 
.3*>9743 + Oo 



- 

1 • 

vP I ''G+GO 

V p 1 a - , n 7 “ G ♦ 3 D 

VP Za 
vn 2a 

TTOT^fftB + Tl 
.2761 VH+31 

"VP3'« 

VP3* 

• 2 7 1 8 h ♦ 3 1 

• 7 90 ''6 1 ♦ 03 

AREAa 
ARE A- 

“^,399743 + 00 
.399743+30 


5^ S 

1 f. . 

1 7 . 

V 1 “ - p 8 7 ^ > r 0 

V p 1 * - - M 7 G : ' ^ : 0 

VP 2* 
VP2» 

“'.276199 + OT 
.232 1 b 8 + 0 1 

VP3- 
vP j* 

- • 7 9 C 6 1 ♦ 0 
-.202188+31 

ARLAa' 

A«t A- 

,399743+30 

.3997H3+0C 



IP. 

19. 

VP|- “.P7^.'"V..C0''^ 
vPi» -.n7b;G3^^C 

VP2 = 
VP?* 

'• 7 8‘B:6 r+03 
- . 7 M ^ 6 1 ♦ 0 3 

■ vP3* 
VP3* 

-.276194+31 

-.276194+31 

are A«“ 

ARE A- 

.3^9743+33 

.399743+3C 



'"2C. " 

2 1 • 

^ P 1 ^ ^ 7 G 0 ‘ ♦ G 

vP I a - , M / " j + GO 

VP 7 a 
V--2* 

^.207l8A+OI 

-.276l9M^0l 

VP3a 
VP Ja 

-•202188-31 
-.74 0 ''61 + 33 

■ are A- 

ARE A ■ 

■,3777‘i3*3c ~ 
,3997H3*')ri 



22. 

V »■' 1 ■ ’ - . 7 ^ ‘ ; 0 

VF’ ?• 


VP Ja 

.740 :6 .♦'^3 

' A RE A - 

,3'»V7H3»7n — 
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I 

o- 


23. 
28 , 

2^. 
2 6 • 
27 . 

VP 1 » ■*>. 87bon3^CD 

VP z« 

-.20218P+01 

VP3* 

*202 1 da+01 

ypi* "•n7‘D0"'w‘*'OC 
vn* -,ii2 5n:+:i 

VP2*» 
VP 2» 

-.78JC61+00 
. 7 |5796**-00 

yP^m 

VP3* 

»27AJ9m*0I 
•?*7 1 39*31 

“■■'vP 1 ■ '- ,1 1 2S'* J + b 1 
vpi- «,ii?5“'C’^rii 

VP 2“ 
VP 2 p 

.|9S559*0I 

.2A7139+01 

VPjn 
yP3- 

« 195S&9+3I 


p i 3 - . 1 I 2 s ^ .i + n 1 

VP 2 * 

•767139+01 

VP3» 

•*.7 1 5796^30 

2 V • 

vP!= 

VP ?*= 

. 1 7S559*0I 

VP J- 

-. 1 95559 + 0l__ 

J f» • 

V P 1 * - . 1 1 7 5 i 

Vp2=» 

,rU797*on 

vP J* 

-.267 1 39*ni 

^ 1 • 

VPI® -«'12S'J^01 

VP 2* 

-..715796 + 0: 

vP3* 

-•?67l39^0l 

J 7 • 

VP 1 «s - , n 75'* ^ + ' 1 

VP 2 » 

19 5559 + 01 

VP3* 

-.1 95559+01 

33 . 

VP 1 a . . 1 1 75 “ ' M 

VP?» 

-.267 1 39+01 

VP3« 

•^.7l 5797 + 00_ 

J . 

v'p 1 i'7 5'-^*'yr 

VP 2» 

^-•267139+01 

VP3* 

*7 1 5796^00 

"i . 

V P 1 ® - . M 7 5 ^ > M 

V p 2 ■ 

19t|Sb9 + 0l 

VP3« 

. 195559 + 0 i 

36 • 

V P i * - . M 7 S !. j + ? 1 ^ 

VP 7* 

- . 7 1 5797 + 00 

VPia 

.2671 39+nl 

3 7 • 

VKI* -.13 75‘w^O ^ 

VP 2» 

. 69 1 53? + 0: _ 

yPjm 

•258o83>OI 

3 H • 

V P 1 « - . 1 3 75' f" 

V p 7« 

• 1HR93';'.CI 

VP3« 

. 18893 r^Ol 

^9 0 

V P j c - • 1 37S'' 1 

VP 7* 

.258 .8 3 + 01 

VP3« 

•69I53?^13 

** U • 

’ 1 i" '-Ti 3 7 5 : -> + 0 1 

V p i * 

• 3^0 1 

VP3- 

-.69|S32*3 3'“ 

8 1 . 

VPl« -.1375:^+01 

VP 7 * 

. 1 fl 4 9 3 r + 0 l 

VP J* 

-•I8fl93:>nl 

•4 7 • 

vP 1 s - - 1 375' 7 4. j 1 

VP 2» 

" ,69|537+0Q^ “ 

VP3- 

-•258ne3+^l 

‘t 3 • 

Vp|* + r>| 

VP7“ 

-.69|532>00 

VP J« 

-.250:83+01 

8 M , 

V P 1 ■ - . 1 3 7 5 ^ J ♦ : 1 

VP 2* 

|qrt93: + G i 

VP3«' 

-• 1 889 3r. + C 1 

M , 

vP|o -.1375^+rl 

VP7 = 

- • 258G8 3^0 1 

VP3« 

-.691 532+00^ 

6 • 

VP 1 a"” - . ! 3 7 5 ' ^ + C 1 

V p 7 » 

— . 75h 0^3 + 0 1 

VP3- 

69I532^00 

/ t 

vP 1 » - . ! 3 7b ' j + 0 1 _ 

V P 2 * 

1 8P93r + (}| 

VPi» 

. 1 88?3i':^0l_ 

Ll • 

V P 1 *' - , 1 3 7 5 '^ ,. + C i 

VP2» 

-.69V5i?^or> 

VP J* 

• 25Br83^or 

8 9 * 

vP 1 3 * - . 1 675 ' l'. + 0 1 

VP?» 

.667?60+0^_ _ 

VP3- 

• ?89-*78+0I 

1. 0 t 

VPl* -,|ft7b'^+''l 

VP 2 = 

.1 fl?30 ! +01"^ 

VPi» 

. IS23CI^0I 

'■ 1 . 

VPI*^ -«l6?5''^+i.:i_ 

V P 2 » 

« 28 9C 7 9 + OJ 

yP3» 

•667268^00_ 


5 2*^ 

VP1» -•!62b'^+:i 

V » » 2 « 

. 2 8 9 7 n + c 1 

VPi* 

— •667260*00 

5 1 • 

V P 1 * -. 1 675"*-+:! 

VP7* 

. 1 <12331 +01 

VP3« 

1 8230 1 *0 » 

5 H . 

V p 1 « - . 1 6 2b: G : 1 

V p7 » 

.667268+00 

VP3» 

- • 28 902H+-P ! 


vpi« -.i6?b':4-:i 

VPZ** 

-.667269+00 

ypjm 

- • 2890784-O 1 

56 • 

‘ \ p 1 ■ - , 1 6 2 5 ^ ♦ c i 

VP 2* 

- • 1 8 2301 +0 1 

VPi- 

1 823C 1 +-31 

b / . 

VP 1 * - . 1 62b'^G^: 1 

VP2» 

-.2890?R^01 

VP3* 

-.66726h + *30 

bH. 

- C a"”^. 1 6 7 b “ ; + c 1 

VP 2a 

-.28902P+01 

VP3« 

• 667268*>Ci0 

57 . 

vpi- -,i67b'-.. + :i 

VP2 « 

187301 >31 

VP3« 

.182301*91 

6 0 • 

VP|* -.167br«.^C.I 

vP 2» 

-.66726P + 0C' 

VP3* 

.289CZfl+0| 


VPl* -♦|Ji7b'J+:i 

VP2a 

. 6830C8+00 

VP3» 

•239972+01 

6 2. 

vpi= -.1075:^+61 

V P 2 » 

,176<>72*0I 

VH3» 

. 175672 + nl 

63. 

VP1» -,IR7b'J+Cl 

V p 7 » 

•239972+01 

VP Jfi 

• 68 3008 + 00 


V P 1 * -Vi 5 75 : j ♦ :T" 

vp2» 

.239^7‘? + 0l 

VP3a 

-•6830C3+00 

65 « 

vPl» -.1875 

vr 2* 

. 1 7 56 7 ?♦ 0 1_ 

VPi* 

•. 175672+01 

6 6 * 

V p 1 ” -^1875:3 + 01 

Vp2* 

♦ 68 3*0 0 8 + 00 

VP3- 

-•239972+01 

6 7. 

ypja -.i87b':G^:i_ 

VP 2“ 

-.683008+00 

VP3« 

-.239972+01 

6 8 . 

" " V p 1 * - . 1 B 7 b r G + 0 1 

VP2* 

"GG I 7 5 677 + 01 

VP3* 

-. 1 75672 + 01 

6 V . 

V rM « • • * 

V p2» 

-.239977+01 

VPJ* 

-.683':08 + 00 

* 7 6 

VfM“ -.107 5:^ + 0 1 

VP7» 

-.239T72*Q 1 

VP3- 

.683303+00 

7 1 ♦ 

V FM » - . 1 8 7 5 0 4. 1 

V P 7 * 

-. 1 75677 + 01 

VP3- 

,175672+01 


V P 1 = - . 1 H 7 1 

V P 2a 

-7683308+00 

VP3a 

.239972+01 

■73. 

V p ! » - . 7 1 2 b - ^ + : 1 

VP2" 

.6 1 8739 + 0:^ 

VP 3* 

.23091 7+01 

7 8 • 

vpi- -.?i7S:^ + ci ' 

VP 2a 

^169383 + Or 

VP3» 

. 1 69''83 + 01 

7b. 

vp|» -.7i?b'': + :t_ 

VP?» 

.23'9 17+01 

VP3a 

.6 1 3739 + 00 

'/ 6 . 

;piG -.2 1 

VP2* 

9 17 + 01 

vP3- 

-.61B739+T0 

7 7 . 

V P 1 = -•7175^^:01 

^ V p 2 ^ 

♦ 1 6 7 : 8 3 ♦ 0 1 

_ ^ vP3* 

-.169083+01 

7 3 . 

‘ V !M » 2 1 7b - b ♦ ‘ 

P 2 a 

.618739.00 

VP 3- 

- . ? 1 5 9 1 ; + 0 1 


• 3^7H3^00 
•J866374.0q_ 

, 38 A637*00 

6 6 J 7 

\3flA6jT7^0C 

• 3866 374 .OQ 

•3B6637^?c . 
•i066i7^?O 
,3S6637^?c;_ 
■; 386&37*CC 

•306637^rtj^ 

'.3866i7^ric 

^373S3I 
'.373b3!>00 
,373S31^Dc_ 
",373b31^?G 
.373S3 1 
•373S31^00 
,373S3 1 ♦OC*^ 
373bJ»*bo 

.373531 
. 373531 ♦bo 

_.37353l^0O_ 

. 37353 1 ♦OC 

_.36nH2M40C^ 
.360‘*2**^bc^ 
__.36«h280C_ 
, 3<J»0‘*2‘^>00 
, 360428 + 00 
36082 8 + 00 
_.36o828+00_ 
.36C828+0C 

,360828 + O 0 

“^r3^b828^0c 
, 36o828+Oo_ 
■" . 36 082 8 + 00 
_.387318^0C_ 
,3873|8^3C 
.38731 

"i’8 ;3i6^oc 
_.3873|8» 0L 
♦38731 8>0C 
•3873|8>0C_ 

^3873 feVco 

. 387318 ^00 
;3873 1 

_ .3873l8^0p_ 
.38731 8>00 
_ .33821 7^0G_ 
.338212^0C~ 

.3382I?^0u^ 

.‘3 38 2 i 2^br 
, 3 Jm2 I 2 + 00 
,33821 2+OC 
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79. 

VPI« -.2I2SCJ*CI 

VP 2 » -•418739400 

VP3» 

-.23091 7*0« 

AR£A« 

•33m2I2"00 

flO’i 

- - VPT» -.71?5?T*or 

VPZ» -•!69C*^3^01 

^VTO’S' 

“im 69C**3*0F 

AREA" 

^X3PTT^^do 

B 1 . 

VPI» 

VP 2 « -•23^9174.01 

VP3- 

-•61 87394-00 

AREA" 

.334212+00 

" B 2 . 

y fj 1 , 2 t‘?5 

V p 2« -.2 J J^TTiTJT 

vPX- 

06 1 07 394-00 

apTa- 

,T3 4 2 I 2^00 

83 . 

V P J » - . 2 1 ? ♦ C 1 

VP 2 - I6908340I 

VP3« 

• 1 69?43^31 

AP£A« 

,33h21 2.00 

Avr 

vPlo -,7|?S^^+Cl 

vp 2 » -•618739+OD 

V^3« 

• 2309 1 74-0 1 

aR£a- 

#3-^4212 + OC 

8 S. 

VP|« -#2375 ■'J, ♦7* 

VP?«" .S9ftft7S + 0 ? 

VP3« 

•221861 * 0 \ 

AREA* 

.321 105.00 

ft 6 • 

““ VP|* -.237S^::*3I 

VP2» .U2M|4^0l 

vPi« 

•162414+01 

area* 

• '3 2 1 1 08 ♦ Ob 

87. 

vPi« -..237s::4.'ri 

VP2» .221 061 *01 

VP3» 

.S9447&+03 

area- 

• 32 11 OS + Otj 

8ft# 

■“ VPt» -.23 7S^u*^3l 

vP 2« •2?I^6I*0I 

VP3« 

-.S94ft7&+00 

AREA- 

• 32l 105 + 00 

y 9 . 

vPi» -#7 178 -:<p.ii 

VP2» #16241 4+0 1 

VP3- 

-•I6241ft+0I 

area- 

.321 IC5 + 0C 

9q; 

VP|» -.737‘:>*‘j*3i 

VP 2» • ft 7S + 00 

VP3« 

-•22IB6|«0I 

ARE A« 

.3?i ios.b'c 

V I . 

V p 1 • - . 2 3 7 s : . ♦ 1 

VP 2* -•b9ftft7S*30 

VP J« 

-•221861+01 

AREA- 

.32! lO^^Oc 

92. 

VP 1 ■ - . 2 3 7S’:'0'^ 7 1 

VP2» -.I6?41ft4-Cl 

VP3» 

-•162414+31 

area- 

.321 !05 + 3o 

93. 

9P|» -.237S'’^ + p| 

VP2« -•221861+31 

VP3« 

-•594479+00 

AREA- 

.321 105 + OC 

9 H . 

VPI* -.23/S^^^rr 

VP2* -*221861+31 

VP3» 

•594475+00 

area- 

.321 lcS+00 

9^ . 

\/P|- *.237S-^ + 0l 

VP2« -.1 624 1 4^0| 

VP3- 

• 1 624 I 4 + 01 

area- 

.321 105+00 

76. 

VPi« -#237b?^-0l 

vp2» -•S9ft4 7~S + 03 

VPjV 

'""•2?lB61 + 0r 

area’- 

7Jn 105 + 00 

ON CONN 

TWaI’IT DI3K NO.i “ 







M r A 

A w£ A 





' 1 

, 1 ' ■ 

. 1 6 3 A ? 5 - 1 





? 

. 3 7 S : C 

•ft9VH74-:i 





3 

7^ "• 7‘r 

• bT "8 ) ? 3 - 5 1 





ft 

.fl7s:?:>0T 

• 1 1 4s 374. 3r 





■ 

■"#'! ! 2^1. 31 

• I 47262 + : 3 





h 

* 1 1 

• 1 797fl7 + :*r 





7 ' 

♦VI 

. 2 1 2 7 1 2+ 





H 

, 1 ft 7 S « ? ♦ '2 1 

.24^4 3 7 + -An 






.7r2'>'‘^»Ti 

. 27d I 6?* :o 






OU I 'iT 

are:a 


1 ■ " .3h:>m 2* 3c 

2 * 376337^}0 
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viEiw Factor output 


ruT AL sA ffP L r’g jTTKTTf NS T A Rf « "T 2T *r ^ fTHTTF”. 0 ^nrrnr“* TTS^ 

EAKTA « I ~ — 


JKkC.tl M(i« I 


rA^vr,F:T main surface, L • I 


U1 

t 

LO 

00 


"A;?n~ 

170 . nrr- 

area 

mrw 

MITAREA — 

AKTA 

91 rWT ACTOR BY 

N 1' M M F R 

H I TS 

NUMbER 

factor 


NUMBFR 

RECIPROCITT 

1 

) 3. 

I 

M. 3333-01 

Z.S32fl-'Cl 

1 

8.S998-0H 

. 

qM ” 


— \ . 

~77532fi-ffr 

2 

2.9107-03 

3 

1 9 • 

3 

1.HA67-33 

2.S378-0 1 

3 

6.8798-03 

— jjj — 

■*“ 1 

M 

iv^7:^=^Tr" 

7 ■ bJ2^-0 

M 

6. 7475—03 

S 

m9 • 

S 

1 .6333-03 

2.S328-0I 

s 

3.2414-03 


■ flV — 

^ & 

■;7.<!>6‘67-yi 

2.53?d-C 1 

6 

5.2922-01 

7 

0. 

7 

**• • ; . ' i 

2.6328-C 1 

7 

o.cocc 

a 

c; 

d 


'2 *'S3 2d-C r 

8^"' 

o.cooo 

9 


9 

*>.L '’.?r 

2.S328-0I 

9 

0 . coac 

0 

c* 

P’i 

. V •/ 

— ?Tsm-arr 

ro' 

0 . COOG 

1 1 

0. 

1 1 


2.S328-C 1 

1 1 

o.ooco 

12 

?. 

12 

’~oT3rT? 

2.S328-C 1 

12 

' o.cooo 

1 3 

1 3, 

1 3 

9.3333-04 

2.02lA-*:i 

13 

1 .0771-03 



— rr. 

j-q 

7 HTb ^TT -^^crr" 

— ?“.D7r6-rr 

■ m 

6.3818-C*3 

\ ^ 


IS 

1 A67-03 

2.02 1 6-Ci 1 

IS 

1.2846-02 

1 6 

f7A. 

F5 

S.tiAA7-53 

2.02 \ *-0 1 

16" 


1 7 

70. 

47 

7,3333-03 

2*02I6^\jI 

17 

5.8016-03 

\P^ 

9',^ 

i’d 


“ 2'.02 l i.-c r 

18 

7.4592-04 

1 9 

C . 

19 

^ " 0 

2.0216-CI 

19 

O.COOO 



n — 

2"C 


— rro2iT-in 

20 

O.COOO 

2\ 

?• 

2 1 

'3*c: 3? 

2.0216-C » 

21 

0.0300 

2 2 

— 

22 


7.52 16-Cr 1 

22 

o.cooo 

2 3 

0 • 

23 

r . 0 '': M 

2.02IA-dl 

23 

0.0000 

- 

" f.” 

2'M 

. j 7 ? n " ■ 

— 2 , 02 l 6 -cr 

29 

“ 0.0003 


2S 

9. 

2S 

3 . "'r.-OH 

1 .5ir»4-Cl 

25 

9,9839-04 

'2<l7 ' 


26 

. 

.>1 

tl 

~fV5ip4-c i 

26 “ 

1.1981-02 

2 7 

1 94. 

' 27 

6. 4A67-03 

1 .5104-C! 

27 

2. I52I-C2 

? n 

2l 3. “ 

2b 

7 . 1 -3C'"-"3 ■ 

1.5404-01' 

28 

3.3621-C2 

29 

1 r ^ . 

2 9 

3.5331-33 

1 .5404-0 1 

29 

1.1759-02 

.. . - ^ 

1 A . 

■3 5 

^r.”3 iT3 - 5'4“ 

1 .5 1 o4-c r 

30 

1 .7749-03 

3 I 

, • 

3 ) 


1.51 r 4-c 1 

31 

O.OOMri 

^ ? 

1 • 

JZ 


1 . 5 1 04 -Ci 1 

32 

■" 5.0 3^3 
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“Ir 


:l'i 





in" 


0* 

33 

n.croo 

1 .SIOH-OI 

33 

o.cooo 


34 

r • 09 Or, 

I.SlCH-Ql 

34 

0,0000 

3S p. 

3S 

• C "* ? D 

I.SIPR-CI 

3S 

0,0000 


“ ' 3 & 

^.C'^on 

1 • 5 1 0 4 - C I 

36 ‘ ' 

“ 0,0000 


total h t r • 

rH^6 • 

Factor 

« 4.8867-02 

RECIPR 

vr - 3,^7SI-03 


T A ^ L r 'JO, i 


ON CHNSTKAtNT DISK PASSING THRU Pl* L ■ 2 


Ul 

I 

00 

vO 


Q 

'5'-€ 

I>. c3 

O 

fc l 


o 

Hi! 


A K A 

NO. or 

ARt A 

View 

HI tawea 

AftEA 

VlEiA factor by 


y R 

HITS 

N U M H E R 

factor 


number 

reciprocity 


. I 

0. 

1 

^ ' n 

b,30l 4-02 

1 

0.0000 


? 

? , 

2 

• G ' 0 

S.3CI 4-02 

2 

0,0000 


3 

?• 

3 

.'0 

b,30I 4-02 

3 

0,0000 


4 

0 . 

4 


5.331 H-C2 

4 

o,cooo 


s 

?. 

b 


5 . 3? 1 4-0? 

S 

o.cooo 


6 

0. 

6 

: '7 n 

b, 331 4-02 

6 

o.rooo 


7 

9. 

7 

■? . C " JO 

b# 30 1 4-C? 

7 

o.aoo? 


a 

“* 6T ■ 

a’’ 

*’.0 

b* 30 1 4-C2 ' 

"a 

0.0000 


9 

r . 

9 

^ . 'J *' '*1 n 

b.3': 1 4-02 

9 

O.fHOO 


1 6 

?• 

iO 

?p 

b.33 1 4-u2 

id 

P,C oco 


1 1 

c. 

i 1 

- c : 9 

b, 30 I 4-G2 

1 1 

o.cooo 


12 ‘ 

37' 

Y2 ' 


5. 36i H-az 

12 

0.0000 




hTT ■ 


vifa' Factor 

* 0,0300 


RECIPR VF - o.cooo 



I A R fi t T NO. I 


Ofj constraint disk Passing thru p?* l * 


3 


area 

N U 4 M F R 

'NO. or ■ 

MI TS 

Af<t A 
NUHdER 

vrEis 

factor 

hTTa^a 

~KW€k 

number 

VlT. f'actoR by 
REC I PROC I TY 

, 

?• 

1 


2.6iaC"Ci 

I 

OfOOOO 

? 

Of 

2 

*> • C '*■ .? 7 • 

2,6100-01 

2 

o.cooo 

3 

0. 

3 

A • ... ‘ ■ ^ t 

2,6 1 aj-G 1 

3 

0,0003 



i ' 

-1. 

2.6 1 8G-3 ! ' 

‘4 

c.oooo 

s 

n. 

b 


2.6I0L--C1 

b 

C.C030 

"^6 

Of 

“ 6" 


7.6l0.i-C 1 

6 

0.O03C 

7 


7 

^ . r 

7 , 6 1 H - C 1 

7 

0 . 0 3 C ? 

8 

9 • * 


- . -r *• : r 

, 6 1 a J - g 1 

0 

0.C3Q3 ' 


II iiir 
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9 

0* 

9 


?.6l«C-0! 

9 

0*0900 

ro * 


^ 


?.6l8C,-C 1 

10 

0.0000 

1 1 


1 1 


2.AI0O-C! 

1 1 

0.0000 

1 2 

b* 

(2 

n • or or 

T76 fad-c 1 

12 

0.000*0 


total rt r*. — n-; vtew rATTOT”- REcTPR'~^r-- ■ o;t5fJO'*r 


Ui 

I 

4^ 

O 
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T A RGCT 'iJO ; z 


Ta^gTt Ta In' SUR^AdTi L « i 


A-it A 
u M tj r R 

UO* OF 

A «L A 

V I EiA 

HITAREA 

area 

VlE«l FACTOR BY 

M I T S 

NUMUE W 

factor 


number 

fiEC ! PROC I f Y 

\ 

S3, 

\ 

1 .M 33 3-D3 

4. 1 28S-0 ! 

1 

l.TRSl-PT 

2 

. 

2 

!.fl667-j3 

4* I2RS-CI 

2 

2.2727-03 

3 

R 1 . 

3 

2.7^ :n-D3 

4, I2flb-ul 

3 

3.2B73-P3 

H 

1 : 1 . 

*+ 

3.3ft67-J3 

4, 1 205-C 1 

4 

4.C"?9n-C3 


• 


2,2667-:J 

4 . 1 PflS-C 1 

S 

2.7S97-03 

h 

?9, 

& 

9. 66i>7- 3*^ 

4 • 1 201^-0 1 

6 

1 * 1769-03 

7 

13. 

7 

H. 333 3 -34 

4,I20S-J 1 

7 

'5.27&9-.04 


1 . 


3.3333-3S 

4. 1 24S-C1 

S 

4 •OSB4-OS 

'" 9 ' ’ 

0 • " 

.. 


4* I28b-G1 

9 

0.C3C0 

10 

c. 


*' • r ' r 

4, 1 20b-C 1 

10 

o.coco 

‘1 ! 

1 . 

il 

3.3333 -Ts 

4 . I 20S-C 1 

1 1 

4.CS94-05 

1 2 

!7. 

1? 

c; , 6667-34 

4. 120S-CI 

12 

6.B993-04 


1 3 

42 . 

1 3 

1 .B' -03 

3.9974-0 1 

13 

1 .7604-03 

\ 4 

0 r • 

i4 

2.7* :--:3 

3.997H-0 » 

14 

3.39&1-C3 

! 

1 05. 

15 

3 . 5 : 3 r - 0 3 

3.9974-01 

IS 

4.4DU-03 

1 6 

1 1 . 

16 

3. 3667-33 

3 . 99 7 H-C 1 

16 

4.2334-C3 

1 7 

06 • 

1 7 

2.0667-33 

3.9974-0 1 

17 

3.6047-03 

1 0 

3 3. 

1 H 

! . 1 ''^..-03 

3^.9974-0! 

18 

1 .3032-03 



' '1 9 

1 .6667-34 

3.9974-0 T 

~ 1 9 

2.C9S7-04 

20 

0. 

2? 


3.9974-0 1 

20 

o.cooc 

2 \ 

0. 

2 r 

^ . c ; : " 

3.9974-G I 

21 

0.0000 

27 

0. 

22 


3.9974-0 1 

22 

0 • C 0 0 3 

2 3 

0 . 

23 

0 . u " 

3.9V74-or 

23 

o.cooo 

24 

4. 

24 

1 .3333-34 

3.9974-01 

24 

1 .6766-04 


25 

?6 • 

25 

0.6667-34 

3.8664-C! 

25 

1 . 1 267-03 

26 

6 0. 

26 

7.2667-33 

3.0664-C 1 

26 

2,9‘)6R-D3 

'~^2r' 

4? . 

27 ' 

\ .4“- :T:- 5^ 

TTS 6 6 4 - c \ 

27 ” ■ 

f .b?aj-03 

2 0 

54. 

2H 

1 .0i 3--T3 

3 .0664-0 « 

28 

2.34CI-03 

2V 

7 1 • 

29 

2.3667-33 

3.3664-01 

29 

3.0764-03 

30 

34. 

' 30 

1 . 1 333-03 

3.0664-01 

30 

1 ,<47JH-03 

3 1 

2. 

31 

6 . 6667-05 

3.0664-Cl 

31 

8,6671-05 

37 

C . 

32 


3.0664-C 1 

32 

o.eoon 

“33 

0. 

33 

- . - ■ )p 

3.8664-0 1 

33 

O.COCJ 


?. 

34 

■2* , 7 or 

3.066 4-0 1 

34 

0.0300 

JO 


3S' 

. A. " 

3.8664-0 1 ' 

35 

O.OQCO 

J6 

2. 

' 

6.6667-05 

3.0664-01 

36 

0.6671-05 


37 “ 

n 

■ ■ 

. 1 . 3333-35 

r;7 j 5 j-c 1 

37 

4.4 0 56-CTS 

J e 

3 . 

30 

! - J4 

3 . 7353-w 1 

30 

1 . 3457-04 

JV 

1 . 

J9 

3 - 3 33 5-^5 

3. 735 J-C 1 

39 

4. 4 056-05 
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*^0 

2« 

HO 

*.66A7-D5 

3.7353-ai 

HO 

a. 9712-06 

Ml 

h • 



3. 7563-01 

4 1 

2.69IH-04 

H2 

2. 

4? 

A«6A67“0S 

3.7353-01 

4 2 

8.9712-06 


0. 

43 

n . C ” C‘ 

~5T7y6'3-Cl 

43 

O.OOCO 

M 

C. 

4 H 


3.7353-CI 

44 

c.oooo 


?. 

MS 

n . ' r ■> 

3. 7353-CTi 

4S 

o.c 000 


c* 

46 

• c r * f . 

3.7353-C I 

46 

O.CODO 

' 'h7 


H 7 ” 

f .7 

3.73S3-CI " 

47 

“ o.roGo 



4 R 

n.c' in 

3 . 7353-01 

40 

0 .0000 


’fV n T9 yTTTT2‘-*5T <1 T 9 . 2 975 


‘->0 

14 6. 

be 

4.0667-03 

3.6'^42-CI 

SO 

6.7871 -C3 


’■■43 V. ^ 

■sr 

1 . HSOr “02 

3.6042-C 1 

5 1 

2.0222-02 

s? 

370. 

S2 

1 . 2333-1? 

3.6042-C I 

S2 

1 .7?0?-t>2 

S3 ' 

rs&". 

' ■ 5 3^‘ 

S.Z-' 'U.-03 

3. <5042-6 1 

53 ^ " 

7.2520-03 

S4 

1 • 

&4 

3,3333-05 

3.6042-C 1 

54 

4.6487-05 

*:.s 

r' p 

bS 

“ • !- • '•r 

“ 3 .6T 4 2-01 

55 

O.cooo 

b6 

0 . 

b6 

' :r 

3.6042-0 1 

56 

0.C30C 


0 . 

-^7 

^ 3 ;; 

3.604^-C I 

“57 

o.nooo 

^0 

0 . 

bb 


3.6?42-C \ 

58 

O.nocc 

S9 

“ ■ Oa' 

S 9 ' “ 

. g: 30 

3.6c 4 2 -Cl 

’59' 

o.cooo “ 

6 0 

0. 

6 : 

4 . c, ' **0 

3.6042-CI 

60 

o.roGo 


<>l ■ T. ffl - o.'irTr y."T73 2“0i ■ ' o.ooco' 


6? 

174, 

62 

5.8' :r,-D3 

3.4732-Cl 

62 

0 . 3943-03 

63 

5u7. 

o3 

1 .623.3-02 

S.'17 32-C 1 

63 ■ 

2.6388-02 

6 4 

52 1 . 

64 

1 . 7367-''2 

3.4732-01 

64 

2.51 34-02 

65 

• 

65 

6. 667 — 03 

3.‘l732-£:i 

6 5 

8.780C-CT3 

66 

0 . 

66 

« ^ r T f* 

3.H732-0 t 

66 

0. COOO 

67 

” C. 

67 

r .0^ 

■~0.4^732-C I 

67 

p , fj n 3 0 

6B 

?. 

68 

•'. 

3.4732-Cl 

60 

O.OOCO 

69 

c. 

(j9 


3 . *1/32-01 

“69 

0.0090 

70 

8 . 

7C 

0 . c ^ 0 

3.9732-C I 

70 

o.cono 

71 

c. 

71 

r*. c 

“T;y7 3 2 -^1 

7 1 

o.ooco 

72 

0. 

72 

n . 0 :• 0 r 
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6. CONCLUSIONS AND RECOMMENDATIONS 

The SRB plume thermal radiation code is capable of defining the 
influence of plume flowfield structure on the magnitude and distribution of 
thermal radiation leaving the plume. Radiant heating rates may be calcu- 
lated for a single SRB plume or for the dual SRB plumes astride the Space 
Shuttle, The plumes may be gimbaled in the yaw and pitch planes. Space 
Shuttle surface geometries are simulated with combinations of quadric sur- 
faces. Surface shading effects are accounted for. The code also has the 
capability to calculate view factors between the SRB plumes and Space Shuttle 
surfaces as well as surface-to-surface view factors. 

The capability of the SRB plume thermal radiation code can easily be 
extended to more accurately model the physics of the radiation transport 
phenomenon. The following e.xtensions in the codes capability are suggested; 

• Model the wavelength dependence of radiative properties on the 
magnitude and distribution of spectral thermal radiation leaving 
the plume. 

• Provide the capability to calculate radiative interchange factors 
between surfaces. This will include the ability to model absorp- 
tion and reflection at the surface. Reflection can be modeled as 
specular, diffuse and real surface reflection. 

• Provide the capability to perform a complete calculation of radia- 
tive interchange factors between all surfaces of a Space Shuttle 
enclosure such as the volume contained within the SRB shroud. 
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